Semi-Inclusive Charged-Pion Electroproduction off Protons and Deuterons: Cross 
Sections, Ratios and Access to the Quark-Parton Model at Low Energies 
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A large set of cross sections for semi-inclusive electroproduction of charged pions (tt ) from both 
proton and deuteron targets was measured. The data are in the deep-inelastic scattering region 
with invariant mass squared W 2 > 4 GeV 2 (up to ~ 7 GeV 2 ) and range in four-momentum transfer 
squared 2 < Q 2 < 4 (GeV/c) 2 , and cover a range in the Bjorken scaling variable 0.2 < x < 0.6. The 
fractional energy of the pions spans a range 0.3 < z < 1, with small transverse momenta with respect 
to the virtual-photon direction, P 2 < 0.2 (GeV/c) 2 . The invariant mass that goes undetected, M x 
or W' , is in the nucleon resonance region, W' < 2 GeV. The new data conclusively show the onset 
of quark-hadron duality in this process, and the relation of this phenomenon to the high-energy 
factorization ansatz of electron-quark scattering and subsequent quark — > pion production mecha- 
nisms. The x, z and P 2 dependences of several ratios (the ratios of favored-unfavored fragmentation 
functions, charged pion ratios, deuteron- hydrogen and aluminum-deuteron ratios for ty + and tt~) 
have been studied. The ratios are found to be in good agreement with expectations based upon a 
high-energy quark-parton model description. We find the azimuthal dependences to be small, as 
compared to exclusive pion electroproduction, and consistent with theoretical expectations based on 
tree-level factorization in terms of transverse-momentum-dependent parton distribution and frag- 
mentation functions. In the context of a simple model, the initial transverse momenta of d quarks 
are found to be slightly smaller than for it quarks, while the transverse momentum width of the 
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favored fragmentation function is about the same as for the unfavored one, and both fragmentation 
widths are larger than the quark widths. 



PACS numbers: 13.60.Le, 13.87.Fh 

I. INTRODUCTION 

There has been a growing realization that understand- 
ing of the resonance region in inelastic scattering, and 
the interplay between resonance behavior and a high- 
energy scaling phenomenon in particular, represents a 
critical gap that must be filled if one is to fathom fully 
the nature of the quark-hadron transition in Quantum 
ChromoDynamics (QCD). The decade or so preceding 
the development of QCD saw tremendous effort devoted 
to describing hadronic interactions in terms of iS-matrix 
theory and self-consistency relations. One of the pro- 
found discoveries of that era was the remarkable rela- 
tionship between low-energy hadronic cross sections and 
their high-energy behavior, in which the former on aver- 
age appears to mimic certain features of the latter. 

At low energies, one expects the hadronic scattering 
amplitude to be dominated by just a few resonance poles. 
As the energy increases, the density of resonances in each 
partial wave, as well as the number of partial waves it- 
self, grows, making it harder to identify contributions 
from individual resonances, and more useful to describe 
the scattering amplitude in terms of a sum of t-channel 
Regge poles and cuts. Progress towards synthesizing the 
two descriptions came with the development of finite- 
energy sum rules, relating dispersion integrals over the 
amplitudes at low energies to high-energy parameters. 

The observation [l| of such a nontrivial relationship be- 
tween inclusive electron-nucleon scattering cross sections 
at low energy, in the region dominated by the nuclcon res- 
onances, and that in the deep- inelastic scattering (DIS) 
regime at high energy, similarly predated QCD. Initial in- 
terpretations of this duality naturally used the theoretical 
tools available at the time, finite-energy sum rules or con- 
sistency relations between hadronic amplitudes inspired 
by the developments in Regge theory that occurred in the 
1960s dS]. With the advent of QCD, De Riijula, Georgi, 
and Politzer offered a qualitative explanation of Bloom- 
Gilman duality Q in terms of either small or on-average- 
canceling higher-twist contributions, but a quantitative 
understanding of the origin of the duality phenomenon 
in QCD remains elusive, although some insight has been 
obtained through phenomcnological model calculations 

Even if it remains counterintuitive that there should 
be a strong relationship between the resonance region, 
in which the lepton scatters from a target hadron tradi- 
tionally treated as a bound system of massive constituent 
quarks, and the deep-inelastic region, where the lepton 
essentially scatters from a single free quark, it is essential 
to provide precise data on the onset of this phenomenon. 
In regions where single-quark scattering is well estab- 
lished, a rich plethora of nucleon structure information 



can be gathered from such reactions through the quark- 
parton model d-Q. 

In this article, we will concentrate on the largely unex- 
plored low-energy domain of semi-inclusive electron scat- 
tering, eA — > e^X , in which a charged pion n^ 1 is de- 
tected in the final state in coincidence with the scattered 
electron. The process of semi-inclusive deep-inelastic 
scattering (SIDIS) has been shown to factorize Q, in the 
high energy limit, into lepton-quark scattering followed 
by quark hadronization. Our focus will be on the process 
where a quark fragments into a pion, such that the elec- 
troproduced pion carries away a large fraction, but not 
all, of the exchanged virtual photon's energy. 

The quark-hadron duality phenomenon has been 
predicted [l0l - iT2l | . and subsequently verified [l3| for 
high-energy meson electroproduction. The relation of 
the duality phenomenon with the onset of factoriza- 
tion in electron-quark scattering and subsequent quark 
hadronization was also postulated and shown to hold in 
e.g. the SU(6) quark model [14E El, El- If so, one 
obtains access to the virtue of semi-inclusive meson pro- 
duction, that lies in the ability to identify, in a partonic 
basis, individual quark species in the nucleon by tag- 
ging specific mesons in the final state, thereby enabling 
both the flavor and spin of quarks and antiquarks to be 
systematically determined. Ideally, one could even di- 
rectly measure the quark transverse momentum depen- 
dence of the quark distribution functions q(x,k t ) by de- 
tecting all particles produced in the hadronization pro- 
cess of the struck quark. A large set of pion electro- 
production data from both hydrogen and deuterium tar- 
gets has been obtained in experiment E00-108 spanning 
the nucleon resonance region. Cross sections for semi- 
inclusive electroproduction of charged pions from 
both proton and deuteron targets were measured for 
0.2 < x < 0.6, 2 < Q 2 < 4 (GcV/c) 2 , 0.3 < z < 1, and 
P t 2 < 0.2 (GcV/c) 2 . The results from this experiment 
permit a first study of a possible low-energy access to the 
quark-parton model, either directly through cross section 
measurements or indirectly through their ratios, possibly 
lowering the energy threshold to access the quark-parton 
model if higher-twist contributions would fully cancel. 

In Section 2, we will describe in detail the relation 
between kincmatical cuts to separate current and tar- 
get region fragmentation events as optimally as possible, 
quark-hadron duality, and a low-energy onset of a factor- 
ized (or precociously factorized) description in terms of 
electron-quark scattering and subsequent hadronization 
of the struck (current) quark. In Section 3, we will relate 
the findings of earlier low-energy experiments to a quark- 
parton model description, and extensions thereof beyond 
the infinite momentum frame including azimuthal-anglc 
and transverse-momentum dependences. Sections 4, 5 



and 6 will cover the experimental details, the data anal- 
ysis procedures, and the systematic uncertainties, respec- 
tively. Finally, Sections 7 and 8 will describe the experi- 
mental results in terms of dependences of ratios and cross 
sections on various kinematic variables, including some 
nuclear dependences, followed by the conclusions. 



II. TOWARDS A HIGH-ENERGY 
DESCRIPTION OF SEMI-INCLUSIVE PION 
ELECTROPRODUCTION 

A. Semi-Inclusive Deep Inelastic Scattering 

In semi- inclusive deep inelastic scattering (SIDIS), a 
hadron h (in our case a charged pion tt ± ) is detected 
in coincidence with a scattered electron, with a sufficient 
amount of energy and momentum transferred in the scat- 
tering process. Under the latter conditions, the reaction 
can be seen as knockout of a quark and subsequent (in- 
dependent) hadronization. 

Fig. Q] gives a schematic picture of this process, includ- 
ing the kinematics. An electron with four-momentum 
(E, k) scatters from a nucleon with mass M (taken to 
be the proton mass M p at rest), resulting in a scattered 
electron with four- momentum (E' , k'), thereby exchang- 
ing a virtual photon with four-momentum q = (y, q) with 
a quark. A meson with four momentum m = (Eh,Ph) 
is produced, with the residual hadronic system charac- 
terized by an invariant mass W . As usual, the four- 
momentum transfer squared is defined as Q 2 = — q 2 and 
the Bjorken variable as x = Q 2 jlMv. The latter can be 
interpreted as the fraction of the light-cone momentum 
of the target nucleon carried by the struck quark. Fur- 
thermore, z is defined as z = (p ■ m)/(p ■ q). In the target 
rest (lab) frame, this becomes z = E^jv, the fraction 
of the virtual photon energy taken away by the meson. 
In the elastic limit, z = 1, and the meson carries away 
all of the photon's energy. Finally, we define Pt to be 
the transverse momentum of the meson in the virtual 
photon-nucleon system. 

At high values of Q 2 and v, the cross section (at lead- 
ing order in the strong coupling constant a s ) for the re- 
action N(e, e f Tr)X can be written in the following way 
(see Ref. |15j). 
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where i denotes the quark flavor and e, is the quark 
charge, and the fragmentation function D qi ^ f - K {z,Q 2 ) 
gives the probability for a quark to evolve into a pion 7r 
with a fraction z of the quark (or virtual photon) energy, 
z = Ek/v. The first part of this formula expresses that 
the cross section factorizes into the product of the vir- 
tual photon-quark interaction and the subsequent quark 
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FIG. 1: Schematic diagram of meson electroproduction. 



hadronization. A consequence of factorization is that 
the fragmentation function is independent of x, and the 
parton distribution function qi(x,Q 2 ) independent of z. 
Both parton distribution and fragmentation functions, 
however, depend on Q 2 through logarithmic Q 2 evolu- 
tion [l6j]. The second part describes the dependence on 
the transverse momentum P t , assumed to be Gaussian, 
and the general dependence [131 °f the cross section in the 
unpolarized case on the angle <j>, the angle between the 
electron scattering plane and the pion production plane, 
with A and B, reflecting the interference terms cjlt and 
cttt, respectively, being functions of x, Q 2 , z, P t . An im- 
portant variable for the analysis is the missing mass M x , 
which is the invariant mass of the undetected residual 
system. Here we refer to this quantity as W (M x ) to 
highlight the fact that it could play a role analogous to 
W for duality in the inclusive case [ll| . If we neglect the 
pion mass, W' 2 is given by 



W' =W 2 - 2zv(M + v- \q\ cos6 qm ) 



(2) 



where v = E — E' and 9 qm is the lab angle between the 
virtual photon momentum \q\ and the outgoing meson 
momentum \p\ . As in the usual inclusive scattering case, 
the square of the (inclusive) invariant mass W is given 
by 



W = M 



(--1) 



(3) 



If we further limit the outgoing meson to be collinear with 
the virtual photon momentum and require that Q 2 jv 2 <C 



(1) 1, we can express W' 2 in terms of z, x, and Q 2 



as 



W' 2 = M 2 + Q 2 (1 - z) (- - 1) . 

x 



(4) 



The quantitative differences between Eqs. (0) and (UJ) arc 
small for the described experimental results, and not vis- 
ible on any of the figures in the remainder of this article. 

For the remainder of this article, we will equate the 
"nucleon resonance region" to the condition that W < 2 
GeV, even if the invariant mass W will be beyond the 
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usually denned resonance region, W > 2 GeV. As can be 
easily read from Eq. @, the larger z the fewer hadronic 
states will be involved in the semi-inclusive pion elec- 
troproduction process, with 2 = 1 the (deep) exclusive 
limit. 



B. Factorization 

If one neglects the dependence of the cross section on 
the pion transverse momentum P t and the angle <fr, the 
SIDIS cross section as given in Eq. ([T} can be written as 



^2qi{x,Q 2 ) D q ^(z,Q 2 ). 



(5) 



(At higher orders one has to worry about gluon fragmen- 
tation functions, but this can be neglected for the e nerg y 
and momentum transfers under consideration here |18f ) . 
The question is how well this factorization into indepen- 
dent functions of x and z is fulfilled in practice. 

Initial investigations of the hadronization process were 
made in electron-positron annihilation and in deep in- 
clastic scattering. By now a wealth of data has been ac- 
cumulated to parameterize the fragmentation functions 
as function of z and Q 2 . It is well known that for the case 
of SIDIS one has to worry about separating pions directly 
produced by the struck quark (termed "current fragmen- 
tation" ) from those originating from the spectator quark 
system ("target fragmentation"). This has been histori- 
cally done for high-energy SIDIS by using separation in 
rapidity, 77, with the latter defined in terms of the pro- 
duced pion energy and the longitudinal component of the 
momentum (along the q direction), 



n 



2 \E r +p* 



(6) 



Early data from CERN [19|, |2C§ suggest that a difference 
in rapidities, A77, between pions produced in the current 
and target fragmentation regions ("rapidity gap") of at 
least A77 « 2 is needed to kinematically separate the two 
regions. 

It has been argued that such kinematic separation is 
even possible at lower energies, or low W 2 , if one con- 
siders only clcctroproduced pions with large elasticity z, 
i.e., with energies close to the maximum energy trans- 
fer [53, H3| ■ Figure [2] shows a plot of rapidity versus z 
for W = 2.5 GeV. At W = 2.5 GeV, a rapidity gap of 
An > 2 would be obtained with z > 0.4 for pion electro- 
production. For larger W, such a rapidity gap could al- 
ready be attained at a lower value of z (see Ref. [21TI22I]). 
For instance, one would anticipate a reasonable kinematic 
separation between the current and target fragmentation 
processes for z > 0.2 at W = 5 GeV. The other issue 
is at which energy scales we can make the assumption 
of independence of the hard scattering process from the 
hadronization process. At low energies we would nor- 
mally view the nucleon as a collection of constituent 
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FIG. 2: (Color online) Relation between elasticity z and cen- 
ter of mass rapidity ncM in semi-inclusive electroproduction 
of various hadrons for W = 2.5 GeV, assuming null trans- 
verse momentum. The band for pions reflects the influence of 
transverse momentum. 



quarks, and the factorization ansatz could break down 
due to effects of final state interactions, resonant nucleon 
excitations and higher-twist contributions (23|, even if a 
sufficient rapidity gap would be established. For this ar- 
ticle, we will simply assume that factorization in terms 
of a hard scattering and subsequent hadronization (called 
kinematical factorization in the remainder of this article) 
holds provided kinematical separation between current 
and target fragmentation is possible, and one is beyond 
the nucleon resonance region, W > 2 GeV. 



To give credence to the latter assumption, we observe 
that in the annihilation process e + e~ — > hX, experi- 
mental data HI El beyond z«0.5atVT = 3 GeV 
(W = 1.94 GeV) were historically described in terms of 
fragmentation functions. The region extends to z > 0.2 
for W = 4.8 GeV (W = 2.84 GeV) and to z > 0.1 for 
W = 7.4 GeV (W = 4.14 GeV). For z > 0.3, fragmen- 
tation functions have also been obtained from data [26| 
on ep — > e'-Ti^A" at an incident energy E = 11.5 GeV, 
with 3 < W < 4 GeV. All of these data are beyond the 
(W > 2 GeV) nucleon resonance region as defined above, 
and seem indeed reasonably well understood in terms of 
a simple fragmentation description. 
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FIG. 3: (Color online) The structure function F2 versus x for 
resonance data from Jefferson Lab and SLAC [2!| [3l[, and 
SLAC and NMC deep-inelastic data, at two different values 
of Q 2 . The solid curves show the GRV parameterization [32] 
at Q 2 = 3 (GeV/c) 2 and Q 2 = 1 (GeV/c) 2 . 



C. Quark-Hadron Duality and Precocious 
Factorization 

At energies W < 2 GeV, it is not obvious that the pion 
electroproduction process factorizes in the same man- 
ner as in Eq. (p}. At energies where hadronic phenom- 
ena dominate, the pion electroproduction process may 
rather be described through the excitation of nucleon res- 
onances, N*. and their subsequent decays into mesons 
and lower lying resonances, N'* [23|, |27j- It has been 
argued that a factorization similar to the one at high- 
energy may appear to hold at low e nerg ies due to the 
quark- hadron duality phenomenon 1 1 CW — ll 21 ] . For that phe- 
nomenon to occur, non-trivial cancellations of the angu- 
lar distributions from various decay channels fl2l [27l . [28j 
would be required to produce the fast-forward moving 
pion at the high-energy limit. 

In the early 1970s, Bloom and Gilman made the phe- 
nomenological observation that there exists a duality be- 
tween inelastic electron-proton scattering in the reso- 
nance region and in the Deep-Inelastic Scattering (DIS) 
regime [lj. More detailed studies over the last decade 
have shown that quark-hadron duality is exhibited over a 
broader kinematic ran ge, a nd with greater precision, than 
was previously known |29l l30j . Duality was found to also 
work quite well locally, with various resonance regions 
averaging to DIS scaling expectations to good approxi- 
mation (< 10%), even down to low momentum transfer 
values (Q 2 « 0.5 (GeV/c) 2 . Alternatively, the individual 
resonance scans average to some global curve even down 
to Q 2 w 0.1 (GeV/c) 2 H [p. This global curve then 
coincides with the DIS scaling expectations at larger Q 2 
(or x). This is illustrated in Fig. [3] The observation 



of duality tells us that higher twist terms mostly cancel, 
or are small, even at these low values of Q 2 , when av- 
eraging over a sufficient (but relatively small) amount of 
resonances and the underlying non-rcsonant background 
contributions. This implies that single-quark scattering 
remains the dominant process, even though one visually 
can see the effect of the quark-quark interactions by the 
resonance peak enhancements. The quark-quark inter- 
actions modify the measured spectrum, bound to create 
confined quarks, but do so only locally. 

Duality studies in inclusive scattering have been ex- 
tended to spin-structure functions, which was predicted 
from both perturbative (33j and nonperturbative QCD 
arguments fl2l . [Hj]. The first experiment accessing 
the spin-dependent asymmetries was SLAC experiment 
E143 [35l437j . Spin asymmetry data reported by the 
HERMES Collaboration (H and JLab (CLAS and 
Hall A [4(j) in the nucleon resonance region were also 
found to be in reasonable agreement with those measured 
in the deep-inelastic region [§3, l4l| - |43j , with possible ex- 
ceptions in the N — A resonance region. More recently, 
both CLAS and Hall A Collaborations have accumulated 
a set of precision data to study the onset of quark-hadron 
duality in polarized inclusive electron-nucleon scattering 
as function of Q 2 |44| , with good agreement found at 
larger Q 2 (> 2 (GeV/c) 2 ). 

While the phenomenon of duality in inclusive electron 
scattering is thus well-established, duality in the related 
case of semi-inclusive meson electroproduction was not 
experimentally tested before this experiment. To exper- 
imentally investigate the existence of quark-hadron du- 
ality in semi-inclusive pion electroproduction processes, 
and how this may be related to a precocious (low-energy) 
factorization and partonic description, was one of the 
main goals of the E00-108 experiment. 

Carlson [l(| suggested several phenomena one could 
look for to explore any possible dual behavior between 
electroproduction of mesons in the resonance region and 
the high-energy scaling expectations by using 'meson tag- 
ging' in the final state, in close analogy to the original 
inclusive case findings by Bloom and Gilman [l|: 

• Do we observe scaling behavior as Q 2 increases? 

• Do the resonances tend to fall along the DIS scaling 
curve? 

• Does the ratio of resonant to non-resonant strength 
remain roughly constant with increase of Q 2 , as was 
one of Bloom-Gilman's original observations? 

Existing experimental charged pion electroproduction 
data show hardly any nucleon resonance structure at 
W > 1.4 GeV, and seem to scale [45| . hinting at par- 
tial answers to some of these questions. In addition, an 
initial investigation of duality in semi-inclusive pion pro- 
duction was made in Ref. |12| . where the factorization 
between parton distribution and fragmentation functions 
was found to hold when summing over the N* resonances 
in the SU(6) quark model. 
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The existence of low-energy kinematical factorization 
fl2l |46| in combination with the quark-hadron duality 
phenomenon may very well lead to a precocious descrip- 
tion of the SIDIS process at low energy in terms of the 
quark-parton model. Applying that to the case discussed 
here, one could anticipate that factorization and a quark- 
parton model description work reasonably well for z > 
0.4 and at relatively low W scales (below 2 GeV). 

In this discussion we have neglected the dependence of 
measured pion yields, as in Eq. ([5]), on the pion trans- 
verse momentum, P t . At high energies the dependence 
on P t has historically been described with a Gaussian 
dependence as exp(—bP 2 ), where 6 _1 represents the av- 
erage transverse momentum squared of the struck quark. 
At lower energies, the measured Pt dependence must re- 
flect to some extent the decay angular distributions of 
the electroproduced resonances in regions where these 
resonances dominate. One would therefore expect the Pt 
dependence to vary with W' at low W' . 



III. SEMI-INCLUSIVE DEEP INELASTIC 
SCATTERING AND THE QUARK-PARTON 
MODEL 

In this Section we will first revisit experimental in- 
formation at relatively low energies, in order to see if 
that exhibits characteristics of a factorized description 
as portrayed by Eq. ([5]). Theoretically, such a factorized 
description is only valid at leading order in a s , and af- 
ter integration over the transverse momentum P t and the 
azimuthal angle <p. Then, we will investigate what can 
be learned from a phenomenological description of mea- 
surements of such transverse momenta and azimuthal an- 
gles, where the factorized description breaks down even 
at leading order, and a whole series of further assump- 
tions must be made to relate these data to a quark-parton 
description. 



A. Low-Energy x — z Factorization and the 
Quark-Parton Model 

Several pieces of evidence suggest that factorization 
may hold at low energies in meson-tagged reactions. Ini- 
tially, skepticism existed about the applicability of the 
quark-parton model at energies below those historically 
used at, e.g., the Electron Muon Collaboration experi- 
ment [471 ] . because of the possibility of overlapping cur- 
rent and target fragmentation regions. Interest grew with 
the findings of the HERMES experiment at DESY, where 
an intriguing similarity was found between results from 
semi-inclusive deep inelastic scattering at moderate en- 
ergies 
ergics 



48 1 and a Drcll-Yan experiment at far higher en- 
49j . This similarity suggested that factorization, 



and a quark-parton model description may after all be 
valid at energies where it is not necessarily expected to 
work. 



The HERMES experiment measured semi-inclusive 
pion electroproduction (7* A — > ir X) in the DIS regime, 
over the ranges 13 < v < 19 GeV and 21 < W 2 < 
35 (GeV) 2 , with an average four-momentum transfer 
(Q 2 ) = 2.3 (GeV/c) 2 . The HERMES analysis explic- 
itly assumed factorization in order to extract the sea 
asymmetry d — u. In particular, it was assumed that 
the charged pion yield factorized into quark density 
distributions Qi(x) and fragmentation functions D™, (z): 

N**(x,z) ex ]Te 2 [ qi (x)D^(z) + q^x) (z) 

i 

(7) 

Indeed, agreement was found between the extracted fla- 
vor asymmetry of the nucleoli quarks sea results of HER- 
MES, and the FcrmiLab Drcll-Yan experiment E866 that 
first reported this flavor asymmetry (at dramatically 
higher energies). Revisiting these HERMES data, at an 
average W > 5 GeV, and constrained to fractions of the 
virtual photon energy, z, of larger than 0.2, it is per- 
haps not so surprising that these data support that the 
factorization assumption used in the HERMES analysis 
appears to be valid for the nucleon sea, even at relatively 
low energy loss. A rapidity gap rj > 2, rendering potential 
sufficient separation between the current and target frag- 
mentation regions and thus kinematical factorization, can 
already be attained at values of z of 0.2 (see Ref. (2~il I22I] ) 
for W = 5 GeV. Given that in the HERMES kinematics 
also W remains larger than 2 GeV, both requirements we 
assumed to be needed for a valid high-energy factorized 
description are fulfilled. 

At even lower energies, with kinematics close to the ex- 
periment reported upon here, a series of measurements 
of semi-inclusive pion electroproduction was carried out 
at Cornell in the mid 1970s, with both hydrogen and 
deuterium targets [50Tj52| . These Cornell measurements 
covered a region in Q 2 (1 < Q 2 < 4 (GeV/c) 2 ) and v 
(2.5 < v < 6 GeV). The results of these measurements 
were analyzed in terms of an invariant structure function 
(comparable to JV" (x, z) of Eq. [7]), written in terms of 
the sum of products of parton distribution functions and 
parton fragmentation functions. The authors concluded 
that this invariant structure function shows no Q 2 de- 
pendence, and a weak dependence on W 2 , within their 
region of kinematics, which would be consistent with a 
factorized quark-parton model description. 

This is the more striking if one realizes that the Cornell 
kinematics cover a region in W 2 between 4 and 10 GeV 2 , 
and in z between 0.1 and 1. In fact, if one would calculate 
W', these results are for an appreciable fraction of their 
kinematics in the region 1 < W' < 2 GeV, which is gener- 
ally associated with the nucleon resonance region. Even 
more, the final pion momentum is often only 1 GeV/c, 
such that final pion-nuclcon scattering effects, especially 
differences between both pion charge flavors, cannot be 
neglected. To complete this enumeration, P t , the av- 
erage transverse momentum of the meson, was typically 
less than 0.1 GeV/c. Unfortunately, not enough statistics 
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and information are available to warrant a careful check 
of duality or factorization in the Cornell data, even if the 
data are suggestive that quark-hadron duality in charged 
pion electroproduction, and a precocious low-energy fac- 
torization, may work in these kinematics. 



B. Transverse momenta and azimuthal angles 

A central question in the understanding of nucleon 
structure is the orbital motion of partons. Much is known 
about the light-cone momentum fraction, x, and virtu- 
ally scale, Q 2 , dependences of the up and down quark 
parton distribution functions (PDFs) in the nucleon. In 
contrast, very little is presently known about the depen- 
dence of these functions on the transverse momentum 
k t of the parton. Simply based on the size of the nu- 
cleon in which the quarks are confined, one would ex- 
pect characteristic transverse momenta of order a few 
hundred McV/c, with larger values at small Bjorken x, 
where the sea quarks dominate, and smaller values at 
high x, where all of the quark momentum is longitudinal 
in the limit x = 1. Increasingly precise studies of the 
nucleon spin sum rule j53l - l56l | strongly suggest that the 
net spin carried by quarks and gluons is relatively small, 
and therefore the net orbital angular momentum must 
be significant. This in turn implies significant transverse 
momentum of quarks. Questions that naturally arise in- 
clude: what is the flavor and helicity dependence of the 
transverse motion of quarks and gluons, and can these 
be modeled theoretically and measured experimentally? 

In the E00-108 experiment, we detect only a single 
hadronization product: a charged pion carrying a (large) 
energy fraction z of the available energy. The probabil- 
ity of producing a pion with a transverse momentum P t 
relative to the virtual photon (q) direction is described 
by a convolution of the quark distribution functions 
and pt-dependent fragmentation functions D + (z,p t ) and 
D~(z,p t ), where pt is the transverse momentum of the 
pion relative to the quark direction, kt is the struck quark 
intrinsic transverse momentum, with the condition (57j 
P t = zkt+pt- The "favored" and "unfavored" fragmenta- 
tion functions D + (z,p t ) and D~(z,pt) refer to the cases 
where the produced pion contains the flavor of the struck 
quark or not. "Soft" non-pcrturbative processes are ex- 
pected [57J to generate relatively small values of pt with 
an approximately Gaussian distribution in p t . Hard QCD 
processes are expected to generate large non-Gaussian 
tails for p t > 1 GeV/c, but probably do not play a major 
role in the interpretation of the E00-108 experiment, for 
which the total transverse momentum P t < 0.45 GeV/c. 

Because the average value of 4> m the E00-108 exper- 
iment is correlated with P t , (sec Fig. HJ, we first need 
to study the 4> dependence. The cross sections for each 
target and pion flavor were parameterized in the form 
of Eq. [TJ The assumed Gaussian P t dependence (with 



slopes b for each case) is an effective parameterization 
that seems to describe the data adequately for use in 
making radiative and bin-centering corrections. Small 
values of A and B are expected from non-zero parton mo- 
tion, as described by Cahn [58| and Levelt- Mulders [59| . 
In general, any non-zero parton motion effects, be it kine- 
matic or dynamic, are proportional to Pt for A, and P 2 
for B, respectively |58l-l6ll]. 

The more recent treatment of Ref. [57j similarly gives 
results for A and B that are very close to zero (espe- 
cially for B). Other possible higher twist contributions 
will also be proportional to powers of Pt/^/Q 2 [62l. loll- 
and therefore suppressed at our low average values of Pt . 
Specifically, the twist-2 Boer-Mulders [6J] contribution 
to B is essentially zero in the models of Ref. [HI, l65j . 
For the kinematics of the E00-108 experiment, the value 
of B for 7r + is expected to be positive and could change 
approximately linearly with x, z and Pt from ~0.002 to 
~0.02 GeV/c, see Fig. 6 and Fig. 7 in Ref. 0. For 
7r~, it is expected to be negative and the dependences on 
x, z and P t to be much weaker. In contrast, the values 
of A and B are much larger in exclusive pion production 
than those predicted for SIDIS. 



IV. EXPERIMENT 

The experiment E00-108 [67| ran in the summer of 
2003 in Hall C at Jefferson Lab. An electron beam with 
energy of 5.479 GeV and currents ranging between 20 
and 60 ijlA was provided by the CEBAF accelerator. In- 
cident electrons were scattered from 4-cm-long liquid hy- 
drogen or deuterium targets and detected in the Short 
Orbit Spectrometer (SOS). The SOS central momentum 
remained constant throughout the experiment, with a 
value of 1.702 GeV/c. The electroproduced mesons (pre- 
dominantly pions) were detected in the High Momentum 
Spectrometer (HMS), with momenta ranging from 1.3 to 
4.1 GeV/c. A detailed description of the spectrometers 
and set-up can be found in Ref. (68j . The experiment 
consisted of three parts: i) at a fixed electron kinematics 
of (x, Q 2 ) = (0.32, 2.30 (GeV/c) 2 ), z was varied from 0.3 
to 1 by changing the HMS momentum, with nearly uni- 
form coverage in the pion azimuthal angle, </>, around the 
virtual photon direction, but at a small average Pt of 0.05 
GeV/c (z-scan); ii) for z = 0.55, x was varied from 0.2 to 
0.6, with a corresponding variation in Q 2 from 1.5 to 4.6 
(GeV/c) 2 , by changing the SOS angle, keeping the pion 
centered on the virtual-photon direction (and again av- 
erage P t of 0.05 GeV/c) (x-scan); iii) for (x, Q 2 ) = (0.32, 
2.30 (GeV/c) 2 ), z near 0.55, Pt was scanned from to 
0.4 GeV/c by increasing the HMS angle (with average cf> 
near 180 degrees) (p-scan). The kinematic settings are 
listed in Table H 
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TABLE I: Kinematic settings (z-scan, x-scan and Pt-scan) in experiment E00-108. The electron beam energy E was 5.479 
GeV. The scattered electrons were detected in SOS set at constant momentum 1.702 GeV/c throughout the experiment. 





V 


Q 2 


X 


W 2 


\q\ 




dm 


z 


Pm 


n — n 

w 


deg 


GeV (GeV/c] 


2 


GeV 2 GeV/c 


deg 


deg 




GeV/c 


GeV 2 


28.71 


3 7Q4 




0.32 




A flQ 


11.54 


11.54 


0.37 


1.397 


3.92 


















0.42 


1.606 


3.65 


















0.49 


1.846 


3.35 


















0.56 


2.122 


3.00 


















0.64 


2.439 


2.60 


















n i a 
U. / 4 


z.oUo 




















0.85 


3.222 


1.60 


















0.97 


3.703 


1.00 


25.70 


3.794 


1.85 


0.26 


6.16 


4.03 


10.55 


10.55 


0.55 


2.082 


3.25 


28.71 


3.794 


2.30 


0.32 


5.70 


4.09 


11.54 


11.54 


0.55 


2.082 


3.05 


31.75 


3.794 


2.80 


0.39 


5.20 


4.15 


12.47 


12.47 


0.55 


2.082 


2.82 


34.55 


3.794 


3.30 


0.46 


4.70 


4.21 


13.27 


13.27 


0.55 


2.082 


2.60 


37.17 


3.794 


3.80 


0.53 


4.20 


4.27 


13.95 


13.95 


0.55 


2.082 


2.37 


39.63 


3.794 


4.30 


0.60 


3.70 


4.32 


14.54 


14.54 


0.55 


2.082 


2.15 


28.71 


3.794 


2.30 


0.32 


5.70 


4.09 


11.54 


11.54 


0.55 


2.082 


3.29 














11.54 


13.54 


0.55 


2.082 


3.29 














11.54 


15.54 


0.55 


2.082 


3.29 














11.54 


17.54 


0.55 


2.082 


3.29 














11.54 


19.54 


0.55 


2.082 


3.29 



The 4> distribution as a function of P t is shown for all 
three data sets combined in Fig.SJ Except for the largest 
x-setting in the a;-scan, the virtual-photon-nucleon in- 
variant mass W was always larger than 2.1 GeV (typi- 
cally 2.4 GeV), in the traditional deep-inelastic region for 
inclusive scattering. In order to avoid complications from 
irN final-state interactions the momenta of the outgoing 
pions were kept greater than 2 GeV/c in most cases. All 
measurements were performed for both tt + and tt~ . 



V. DATA ANALYSIS 

The raw data collected by the data acquisition system 
were processed by the standard Hall C analysis engine 
(ENGINE), which decodes the data into physical quan- 
tities on an event by event basis. The main components 
of the data analysis include tracking, event reconstruc- 
tion, determination of and correction for experimental 
and kinematic offsets, particle identification and event 
selection, background estimation and subtraction, correc- 
tion for detector efficiencies, and electronic and computer 
dead times. Many steps of the analysis here are similar 
to the F v data analysis described in Ref. (68[. Below we 
will discuss some of the steps, and will emphasize details 
relevant to the E00-108 experiment. 

Accidentals: Random coincidences occur between 
events from any two beam bursts within the coincidence 
timing gate. The resulting coincidence timing structure 



of random coincidences is peaked every 2 ns (defined 
by the beam microstructure of the CEBAF accelerator). 
The random events under the real coincidence peak can- 
not be identified but their contribution can be estimated. 
The data were corrected for these random contributions 
by selecting a number of random peaks and subtracting 
their average content from the content of the real coin- 
cidence peak. The accidentals were taken for an ^80 
ns interval (for 40 bursts far from the coincidence peak 
(e'ir~ at negative polarity of the HMS spectrometer, and 
e'p and e'7r + peaks at positive polarity)), and the aver- 
age number of accidentals (within 2 ns) was defined as 

N Accidental s 

( 2 ns) = [N Accidentals (80 ns)]/40. 

Electron Identification: Electrons were identified in 
the SOS using a combination of the SOS gas Cerenkov de- 
tector and calorimeter. The gas Cerenkov detector was 
used as a threshold detector with a mean signal of ~7 
photoelectrons per electron. Good electron events were 
selected for a photoelectron (pe) cut N pe > 0.5. This 
cut was chosen to ensure good efficiency over the full ac- 
ceptance, even after accounting for the position depen- 
dence of the pe yield. To determine the efficiency of the 
Cerenkov detector, an electron sample was selected from 
data with the calorimeter cut E ca i/P e - > 0.8, where E ca i 
is a total energy deposited in the calorimeter, and P e - is 
the momentum of the particle in the electron arm (SOS) 
defined by tracking. The Cerenkov detector efficiency is 
then given by the ratio of events with and without the 
Cerenkov detector cut. The efficiency was found to be 
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FIG. 4: Pt distribution of the data from this experiment as a 
function of azimuthal angle <j). 



a cut N pe < 2. 

The separation of pions from protons ( and partly from 
kaons) relies on the HMS aerogel detector [69| . Whether 
or not a particle traversing the aerogel Cerenkov detec- 
tor produces a signal depends on the index of refraction 
of the aerogel material (n re f) and the particle velocity 
(/3 = v/c). The mean number of photoelectrons (for aero- 
gel material with n re y=1.015) was N pe ~7-8 and slightly 
varied with a particle momentum (for z-scan data). The 
aerogel Cerenkov detector efficiency is 99.5 ± 0.02% for a 
threshold cut of N pe > 1. 

Real electron-proton coincidences are eliminated via 
coincidence time cuts in the analysis. In all kinematic 
settings of the experiment the electron-hadron coinci- 
dence time distribution is well described by Gaussian 
with a < 250 ps (in average). In the analysis a cut ±1.2 
ns is used on the e-7r coincidence time. At the high- 
est momentum setting of the HMS (Phms=^-1 GeV/c), 
which is the worst case, there is still about 3 ns sepa- 
ration between electron-proton and electron-pion coinci- 
dence peaks. Even in the absence of proton rejection from 
the aerogel Cerenkov detector, the protons (and partly 
kaons at low momenta) would be removed in the random 
subtraction. 



> 99.8% for a photoelectron cut of N pe > 0.5. 

Electrons deposit their entire energy in the calorime- 
ter peaking E ca i/P e - distribution at 1. Good electron 
events in the calorimeter were selected by applying the 
cut E ca i/P e - > 0.7. This cut removes most of the pions, 
while keeping high electron detection efficiency. The ef- 
ficiency of the calorimeter was determined in a similar 
fashion as for the gas Cerenkov detector. A particle iden- 
tification cut was placed on the gas Cerenkov detector, 
and the calorimeter efficiency was estimated as the ratio 
of events passing the calorimeter cut to the total number 
of events. The corresponding efficiency was estimated to 
be > 99.5 ± 0.1%. The pion rejection factor (the ratio 
of pions with and without cut on energy deposition in 
calorimeter) in this case was ~ 1:20. 

Pion Identification: Pions in the HMS are se- 
lected with aerogel and gas Cerenkov detectors. The gas 
Cerenkov detector's function was to separate electrons 
from negatively-charged pions. The HMS calorimeter did 
not play a significant role and was used only for cross 
checks, and for the gas Cerenkov detector efficiency de- 
termination. 

Pions under the Cerenkov radiation threshold do not 
in principle produce a signal in the detector. However, 
pions may produce (5-electrons, which will result in a pho- 
toelectron number greater than zero. Applying a cut to 
reject electrons may then reject pion events as well. To 
determine the pion efficiency of the Cerenkov detector 
cut positive polarity n + data were used with a calorime- 
ter cut (to take out contribution from positrons). The 
ratio of events passing the Cerenkov detector cut to all 
the events is then the Cerenkov detector pion efficiency. 
The pion Cerenkov detector efficiency is 99.6 ± 0.05% for 



A. Background subtraction and corrections 

Contribution from Target Walls: Events from the 
aluminum walls of the cryogenic target cell were sub- 
tracted by performing empty target runs (dummy runs). 
The dummy data are analyzed in the same way as the 
regular data including the same method of random co- 
incidence subtraction and applying the same analysis 
cuts. The effective charge-normalized yields are then sub- 
tracted from the real data yields taking into account the 
difference in the wall thickness between the target cell 
(0.133 mm) and dummy target (1 mm). In most cases, 
the estimated contribution of the target can to the mea- 
sured yield is quite small, about 2 — 3%. The uncertainty 
in the ratio of the thickness of the dummy relative to the 
target can (2—3%) contributes to a negligible uncertainty 
to the total yield. 

Radiative Corrections: Essentially all of the events 
that "radiate in" to a given bin come from either: (i) 
incoming electrons with a lower actual energy than the 
nominal beam energy, because they have radiated a pho- 
ton; or (ii) scattered electrons with higher energy than 
the one measured in the spectrometer, because they radi- 
ated a photon. In both cases, the value of v at the vertex 
is lower than the reconstructed one, hence z is larger and 
W' is smaller than the nominal value. 

The radiative tails within our semi-inclusive pion elec- 
troproduction data were estimated using the Monte Carlo 
package SIMC. The radiative correction formula coded is 
based on the work of Mo and Tsai [ztj , which originally 
was derived for inclusive electron scattering, but was 
modified for use in coincidence experiments [7l|. Details 
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TABLE II: The values of radiative corrections for x-scan data. 



X 


z 


Q 2 

(GeV/c) 2 


4 

(%) 


(%) 


4 

(%) 


(%) 


0.22 


0.55 


1.59 


4.2±0.4 


5.6±0.6 


4.6±0.5 


5.4±0.5 


0.26 


0.55 


1.88 


4.5±0.5 


5.8±0.6 


4.9±0.5 


5.6±0.6 


0.30 


0.55 


2.17 


4.9±0.5 


6.2±0.6 


5.2±0.5 


5.9±0.6 


0.34 


0.55 


2.46 


5.3±0.5 


6.5±0.7 


5.6±0.6 


6.3±0.6 


0.38 


0.55 


2.75 


5.8±0.6 


6.9±0.7 


6.1±0.6 


6.7±0.7 


0.42 


0.55 


3.04 


6.2±0.6 


7.4±0.7 


6.5±0.7 


7.1±0.7 


0.46 


0.55 


3.32 


6.8±0.7 


7.8±0.8 


7.0±0.7 


7.5±0.8 


0.50 


0.55 


3.61 


7.3±0.7 


8.3±0.8 


7.5±0.8 


8.0±0.8 


0.54 


0.55 


3.90 


7.9±0.8 


8.9±0.9 


8.1±0.8 


8.5±0.9 


0.58 


0.55 


4.19 


8.5±0.9 


9.4±0.9 


8.7±0.9 


9.1±0.9 



TABLE III: The values of radiative corrections for z-scan 
data. 



X 


z 


Q A 

(GeV/c) 2 


(%) 


(%) 


(%) 


k d 
(%) 


0.32 


0.37 


2.31 


1.6±0.2 


3.3±0.3 


2.1±0.2 


3.2±0.3 


0.32 


0.42 


2.31 


2.4±0.2 


4.1±0.4 


2.8±0.3 


3.9±0.4 


0.32 


0.49 


2.31 


3.4±0.3 


5.1±0.5 


3.8±0.4 


4.8±0.5 


0.32 


0.55 


2.31 


4.5±0.5 


6.2±0.6 


4.9±0.5 


5.8±0.6 


0.32 


0.64 


2.31 


5.9±0.6 


7.5±0.8 


6.2±0.6 


7.0±0.7 


0.32 


0.74 


2.31 


7.8±0.8 


9.3±0.9 


8.1±0.8 


8.8±0.9 


0.32 


0.85 


2.31 


10.8±1.1 


11.9±1.2 


11. Oil. 1 


11.5±1.2 


0.32 


0.97 


2.31 


18.3±1.3 


18.5±1.9 


18.3±1.8 


18.3±1.8 



of the implementation are described in Ref. [72|. The 
original formulation of the radiative correction procedure 
used in SIMC was for (e, e'p) reaction. The formula were 
extended to pion electroproduction by D. Gaskell pZ3| . 

As a cross-check, we also estimated radiative correc- 
tions using the code POLRAD. The standard FORTRAN 
code POLRAD-2.0 [74| was written for radiative cor- 
rection (RC) calculations in inclusive and semi-inclusive 
deep-inelastic scattering of polarized leptons by polar- 
ized nucleons and nuclei. The program, which is based 
theoretically on the original approach proposed in the 
Ref. [75| , was created to suit the demands of experiments 
with fixed polarized nuclear targets and at a collider. A 
new version of POLRAD [7(| was created to calculate 
the RC for semi-inclusive (polarized) experiments. In 
this case the cross section depends additionally on the 
variable z. 

The radiative corrections calculated with POLRAD- 
2.0 are in good agreement with SIMC. On average the 
RC's are on the level of ~ 6 — 8% for all our data sets 
at z < 0.7 and reach ~ 15% at z > 0.9. The relative 
values of radiative corrections at our kinematic settings 
are listed in Tables HU [HI] and M\ 

Exclusive Pions: In addition, we subtracted radia- 
tive events coming from the exclusive reactions e + p — > 
e' + 7r + + n and e + n — > e! + 7r~ + p. This required a 
model for the cross section of exclusive pion electropro- 
duction that is valid for a large range of W (from the 
resonance region to W w 2.5 GeV) at relatively large 
Q 2 . The model used in this analysis started with the pa- 



TABLE IV: The values of radiative corrections for Pt-scan 
data. Note, for this set of measurements the scattered electron 
kinematic was fixed at Q 2 =2.31 {GeV/cf 



X 


z 


Pi 














(GeV/c) 2 


(%) 


(%) 


(%) 


(%) 


0.32 


0.55 


0.01 


6.3±0.5 


7.9±0.7 


6.4±0.5 


7.4±0.5 


0.32 


0.55 


0.03 


5.4±0.7 


6.7±0.7 


5.4±0.7 


6.1±0.7 


0.32 


0.55 


0.05 


4.7±0.7 


5.9±0.7 


4.3±0.7 


5.0±0.7 


0.32 


0.55 


0.07 


3.9±0.6 


4.9±0.6 


4.8±0.6 


5.5±0.6 


0.32 


0.55 


0.09 


1.4±0.7 


2.2±0.7 


3.3±0.7 


3.9±0.7 


0.32 


0.55 


0.11 


1.4±0.7 


2.2±0.8 


2.5±0.7 


3.1±0.8 


0.32 


0.55 


0.13 


-1.6±0.9 


-1.0±0.9 


-0.8±0.9 


-0.1±0.9 


0.32 


0.55 


0.15 


-1.4±1.2 


-0.7±1.2 


-1.5±1.2 


-0.8±1.2 


0.32 


0.55 


0.17 


-3.9±1.8 


-3.4±1.9 


-1.1±1.8 


-0.5±1.8 


0.32 


0.55 


0.19 


-8.3±3.2 


-7.9±3.3 


-5.7±3.2 


-5.3±3.2 



rameterization of exclusive 7r + and tt production cross 
section data from [z3 at W 2.2 GeV and Q 2 = 0.7 
and 1.35 (GeV/c) 2 . This parameterization describes the 
more recent data taken at Jefferson Lab as part of the 
Charged Pion Form Factor program [Hi IzE Izl (W = 
1.95 GeV, Q 2 =0.6-1.6 (GeV/c) 2 and W = 2.2 GeV, 
Q 2 =1.6, 2.45 (GeV/c) 2 ) reasonably well. 

While the starting parameterization is appropriate for 
describing exclusive pion production above the resonance 
region, it does rather poorly for values of W significantly 
smaller than 2 GeV. Since no existing model or parame- 
terization describes exclusive pion production both in the 
resonance region and at large W, we chose to adjust our 
starting model by- hand to give good agreement with the 
MAID model [8(J of pion electroproduction in the reso- 
nance region. This by-hand adjustment began with the 
assumption that the longitudinal contribution was well 
described by the starting model, even at relatively low 
W. Discrepancies between the starting fit and the MAID 
calculation were attributed to the transverse cross section 
and were removed by assuming a more modest W depen- 
dence therein. Wc further simplified the model by assum- 
ing that the TT and LT interference terms mostly aver- 
aged to zero over our experimental acceptance so that 
they contributed negligibly to the radiative events. 

We ran SIMC with this modified model for exclusive 
7r + electroproduction on the proton and for 7r + and tt~ 
production on the dcutcron for all our kinematic settings 
(z-scan, x-scan and P t -scan). 

Contributions from exclusive pions were subtracted on 
a bin by bin basis. On average, the contribution from 
the exclusive tail was estimated to be 4-5% for the x-scan 
data, 5-15% for the z-scan data at z < 0.8, and 8-10% 
for the P t -scan results (see Tables El |VH and IVIII) . 

The radiative tail from exclusive events is the domi- 
nant correction for our data at z > 0.8. For z > 0.9 
the contributions from exclusive pions become more than 
50%. 

We also performed an alternative analysis using the 
code HAP RAD The two results agree to within 

±10 — 15% in the relative contribution of the radiative 
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TABLE V: The relative contribution of radiative exclusive tail 
for rr-scan data. 



X 


z 


Q 2 

(GeV/c) 2 


4 

(%) 


k h 
(%) 


4 

(%) 


(%) 


0.22 


0.55 


1.59 


6.1±0.2 




3.7±0.1 


5.1±0.2 


0.26 


0.55 


1.88 


5.2±0.1 




3.5±0.1 


5.1±0.1 


0.30 


0.55 


2.17 


4.6±0.1 




3.4±0.1 


5.3±0.1 


0.34 


0.55 


2.46 


4.6±0.1 




3.3±0.1 


5.1±0.1 


0.38 


0.55 


2.75 


4.2±0.1 




2.9±0.1 


4.8±0.1 


0.42 


0.55 


3.04 


3.8±0.1 




2.7±0.1 


4.9±0.1 


0.46 


0.55 


3.32 


3.7±0.1 




2.6±0.1 


4.2±0.1 


0.50 


0.55 


3.61 


3.1±0.1 




2.3±0.1 


3.6±0.1 


0.54 


0.55 


3.90 


3.2±0.1 




1.9±0.1 


3.1±0.1 


0.58 


0.55 


4.19 


2.5±0.1 




1.5±0.1 


2.5±0.1 



TABLE VI: The relative contribution 
tail for z-scan data. 



of radiative exclusive 



X 


z 


Q' 2 

(GeV/c) 2 


*t 
(%) 


(%) 


(%) 


v D 
(%) 


0.32 


0.33 


2.31 


3.6±0.2 




2.6±0.1 


3.6±0.2 


0.32 


0.38 


2.31 


3.9±0.1 




3.1±0.1 


4.6±0.1 


0.32 


0.44 


2.31 


4.3±0.1 




3.4±0.1 


4.7±0.1 


0.32 


0.50 


2.31 


4.1±0.1 




2.8±0.1 


5.0±0.1 


0.32 


0.55 


2.31 


5.9±0.1 




4.4±0.1 


7.6±0.1 


0.32 


0.61 


2.31 


7.5±0.1 




5.8±0.1 


8.7±0.2 


0.32 


0.66 


2.31 


8.8±0.1 




6.4±0.1 


10.3±0.2 


0.32 


0.72 


2.31 


11.0±0.2 




7.7±0.1 


12.2±0.2 


0.32 


0.78 


2.31 


13.8±0.2 




8.7±0.2 


15.1±0.3 


0.32 


0.83 


2.31 


15.7±0.3 




9.5±0.2 


18.0±0.4 


0.32 


0.89 


2.31 


21.8±0.4 




15.0±0.3 


30.3±0.6 


0.32 


0.94 


2.31 


^90 




^90 


^90 



exclusive tail. Thus, the resulting uncertainty is only at 
the 1% level or less. 

Diffractive p: Some of the detected events may orig- 
inate from the decay of diffractive vector meson produc- 
tion. The underlying physics of this process, which can 
be described as that the virtual photon fluctuates into a 
vector meson, which subsequently can interact with the 



TABLE VII: The relative contribution of radiative exclu- 
sive tail for Pt-scan data. Note, for these measurements 
the value of four-momentum transfer square was kept at 
Q 2 =2.31 (GeV/c) 2 . 



X 


z 
















(GeV/c) 2 


(%) 


(%) 


(%) 


(%) 


0.32 


0.55 


0.01 


5.8±0.1 




4.0±0.1 


3.1±0.1 


0.32 


0.55 


0.03 


6.9±0.1 




4.6±0.1 


3.7±0.1 


0.32 


0.55 


0.05 


7.4±0.1 




5.3±0.1 


7.3±0.2 


0.32 


0.55 


0.07 


8.7±0.2 




5.5±0.1 


9.2±0.3 


0.32 


0.55 


0.09 


9.0±0.2 




6.0±0.2 


8.4±0.3 


0.32 


0.55 


0.11 


9.6±0.3 




6.4±0.3 


9.6±0.5 


0.32 


0.55 


0.13 


11.1±0.5 




6.8±0.4 


10.3±0.6 


0.32 


0.55 


0.15 


10.8±0.6 




9.8±0.8 


11.7±0.9 


0.32 


0.55 


0.17 


16.8±1.6 




11.0±1.8 


15.6±2.2 


0.32 


0.55 


0.19 


15.9±3.3 




22.0±5.5 


22.1±5.5 



nucleon through multiple gluon (Pomeron) exchange, is 
distinctively different from the interaction of a virtual 
photon with a single current quark. Again, we used SIMC 
to evaluate such a diffractive p meson contribution. 

The p(e,e'p°)p cross section calculation was based on 
the PYTHIA [82] generator, adopting similar modifica- 
tions as implemented by the HERMES collaboration to 
describe lower-energy processes [Hj]. Additional modi- 
fications were implemented to improve agreement with 
p° cross section data from CLAS in Hall B at Jefferson 
Lab H. 

The p(e,e'/o°)p cross section can be written as 



>,Q 2 ) = r T (l + ei?) 



M 2 



T yp->pp 



Ml 



(8) 

where Tt is the transverse photon flux factor, R = <jl I <Jt 
is the ratio of longitudinal to transverse cross sections, 

/ A/ 2 \ " 

( M -i+Q'i ) (n = 2 in PYTHIA) is an additional fac- 
tor that accounts for the suppression of the cross section 
from virtual photons, and a lp ~^ pp is the photoproduction 
cross section. The modifications to the PYTHIA model 
implemented for this analysis mimic those implemented 
by the HERMES collaboration: 

1. The calculation of Tt was performed with no high- 
energy approximations 

2. An improved parametrization of R = (Jl/ct 

3. Replacement of the exponent n = 2 with n w 2.6, 
more consistent with lower energy data 



The t dependence of the p° cross section is 
parametrized as 



(9) 



where t' = t — t mm (< for clcctroproduction) and b is 
the slope parameter. Note that at tf = 0, b also impacts 
the overall scale of the forward cross section. The HER- 
MES/PYTHIA model assumed a value of b w 7 GeV" 2 
for all energies. However, CLAS data suggested that this 
constant value of b did not adequately describe the t' de- 
pendence at JLab energies. The model used in SIMC fits 
b as a function of cAt (the vector meson formation time). 
Above cAr = 2 fm, b was taken to be a constant value of 
7.0 GcV~ 2 , while for cAt < 2 fm, b increased from 1.0 
GeV" 2 to 7.0 GeV~ 2 between cAr = 0.4 fm and 2 fm. 

Using the above model, the fraction of events due to pi- 
ons from the decay of produced p mesons was estimated 
to range from a few percent at low z to about 15% at 
z = 0.6, and was subtracted on bin by bin basis. The 
SIMC determination of the exclusive p° contribution to 
the semi-inclusive yield was also checked independently 
using a program and model developed by the CLAS col- 
labration [851 ] . The two calculations were found to agree 
to about the 10% level. 
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Pion decay: Pion decay in flight is included in the 
Monte Carlo simulation. Charged pions predominantly 
decay via — > fj Vu with a branching fraction of 
99.99%. In SIMC the pion can decay at any point along 
its path in the magnet-free regions and at fixed points in 
the magnetic fields of the HMS. The muon momentum 
is calculated in the pion center-of-mass frame, where the 
angular distribution is uniform and the muon momen- 
tum is fixed. The muon is then followed through the 
spectrometer to the detector hut. Like in the experi- 
mental data, the muon is treated as if it were a pion 
in the reconstruction of target variables. In both the 
experimental and simulated data the muons constitute 
a background that is not removed by means of particle 
identification. However, the distribution of the muons 
in various reconstructed quantities is much broader than 
that of the pions. The fraction of pions decaying in flight 
on their way from the target can be calculated from the 
spectrometer central momentum and path length. While 
at low momentum roughly 20% of all pions decay on their 
way to the HMS detector hut, only a quarter of the muons 
fall within the acceptance and pass all cuts. More than 
85% of all simulated muons that survive all cuts origi- 
nate in the field free region behind the HMS dipole. We 
found the muon contamination after applying all cuts to 
be ~ 10% at lowest momentum 1.4 GeV/c, and ~ 2% in 
the momentum range = 3 — 4 GeV/c. The uncertainty 
associated with pion decay is estimated to be <1%. 

Pion Absorption: Some pions are lost due to nu- 
clear interactions in the materials that the particles pass 
through on their way from the target to the HMS detec- 
tor hut. Pions lost in hadronic interactions are largely 
due to absorption and large angle scattering, resulting 
in pions that do not strike all detectors required to form 
a trigger. The transmission of pions through the spec- 
trometer is defined as the fraction of pions that do not 
interact with any of the materials. 

The calculation of the pion transmission through the 
materials is determined by the choice of pion-nucleus 
cross section. In particular, the total cross section, which 
is defined as the sum of all hadronic interactions, repre- 
sents an underestimate of the transmission. This can be 
explained in terms of the contribution of the individual 
pieces to the effective loss of pions. Elastic scattering is 
peaked in the forward direction (small angles), so that 
a large fraction of the elastically scattered pions are ex- 
pected to still produce a valid pion event. In addition, 
inelastic scattering does not necessarily correspond to an 
invalid trigger. On the other hand, a pion that is truly 
"absorbed" will clearly not result in a trigger. Therefore, 
the transmission is calculated from the reaction cross sec- 
tion which includes all hadronic interactions except for 
elastic scattering (a reac = a absorption + & inelastic)- The 
reaction cross section is approximately the average of the 
total and absorption cross sections and the uncertainty 
on the transmission can be estimated from these two lim- 
iting cases. At all kinematic settings of E00-108, the pion 
absorption was estimated to be below 1 — 2 %. 



TABLE VIII: Fraction of K + relative to tt+ (in %) for hy- 
drogen and deuterium targets, for two different cuts on the 
number of photoelectrons (N pe ) in the aerogel detector. 





(K + /tv + ) h 




z 


N pe > 


N pe > 1 


N pe > 


N pe > 1 


0.97 


9.2 


8.9 


11.4 


10.4 


0.85 


5.6 


5.2 


7.2 


6.3 


0.74 


2.7 


2.0 


3.5 


2.5 


0.64 


0.2 


0.1 


0.3 


0.1 



TABLE IX: Kaon contamination (%) in hydrogen and deu- 
terium data at N pe > 1 cut on aerogel. Estimations are based 
on linear fits (Eq. [TO) of the data from Table IVITlI 



Phms 
GeV/c 


z 




K+(D) K-(D) 


3.703 
3.222 
2.803 
2.439 


0.97 
0.85 
0.74 
0.64 


9.3±0.9 0.6±0.2 
7.6±0.5 0.4±0.3 
4.7±0.6 0.3±0.4 
3.3±1.1 -0.1±0.2 


10.9±0.9 0.9±0.2 
8.9±0.5 -0.2±0.3 
5.6±0.7 0.5±0.3 
3.6±1.3 O.liO.l 



Kaon contamination: For low momentum settings 
(Phms < 2.4 GeV/c), which is the case for our x-scan 
and Pj-scan data, kaon contamination is negligible. At 
these momenta the real e-K coincidence peak will be well 
outside the e-7r coincidence peak, so the kaons are elim- 
inated by the coincidence timing cut. In addition the 
HMS aerogel Cerenkov detector can be used for separa- 
tion of pions from kaons. For kinematics with pion mo- 
menta above 2.4 GeV/c, a correction was made to remove 
kaons from the pion sample. To estimate kaon contami- 
nation, the coincidence timing distribution was analyzed 
for z-scan data. It was assumed that the "real" coinci- 
dence peak is a sum of tt+K, therefore this spectrum was 
fitted with a sum of two Gaussian distributions (86j . 

The summary on kaon contamination is presented in 
Tables IVIIII and IIXI The contamination at negative po- 
larity was found to be very small. For positive polarity 
the worst case (about 10% contamination) was found at 
z > 0.85, but a more typical contamination is less than 
2% 

Kaon contamination from the actual data was cor- 
rected for as: 

Y e c x °; r = Y exp x (1-R k ), (10) 

where 

i?f = 0.265 x (z - 0.63) , , 

Rf = 0.294 x (z - 0.63) [ ' 

for hydrogen and deuterium targets. 

Tracking Efficiency and Multiple Tracks: The 

tracking algorithm performs a \ 2 minimization by fitting 
a straight line through both drift chambers. The track- 
ing efficiency is defined as the ratio of events that should 
have passed through the drift chambers and the number 
of events for which a track was found. The fraction of 
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events that should have passed through the drift cham- 
bers is defined by a requirement on hits in a fiducial area 
composed of a particular set of scintillator paddles. The 
efficiency depends on both the drift chamber hit efficiency 
and the tracking algorithm finding a track. 

On average, the typical tracking efficiencies for pions in 
the HMS are better than 95% in most cases and are often 
better when electrons are detected in coincidence in SOS. 
At both positive and negative polarities of the HMS, and 
hydrogen and deuterium targets, during data taking the 
beam current was optimized to keep the HMS rate below 
500—600 kHz. This helps to minimize the well-known 
negative impact of multiple hits on track reconstruction 
in HMS. In many cases the events with more than 15 hits 
in both chambers are caused by projectiles scraping the 
edge of one of the magnets and causing a shower of parti- 
cles, hence multiple hits in the drift chambers. Examina- 
tion of such events and data taken with high rates in HMS 
and/or SOS spectrometers revealed that sometimes the 
tracking code picks as the "best" among multiple track 
candidates a track that is not what a user looking at the 
event display would have picked. To eliminate this prob- 
lem a new code was written that uses selection criteria 
different from the one track criteria [H| . 

With the new code the resulting tracking efficiencies in 
HMS and SOS arc on the level of 97±0.1% and 98±0.4% 
respectively. (With the improved "pruning" algorithm 
we gain about 2-3% useful tracks). The difference be- 
tween HMS and SOS mainly reflects the difference in 
incident count rates. 

Coincidence and Cerenkov detector blocking: 
Another source of event loss is connected with coinci- 
dence and Cerenkov detector blocking effects. The coin- 
cidence time is determined by a clock that starts when 
an HMS signal arrives and stops when the SOS signal 
arrives. Two effects can cause the coincidence timing for 
good events to fail. In the first case, a random SOS single 
arriving before the coincident particle can stop the clock 
too early, effectively blocking the coincidence. A cut on 
the coincidence time will largely remove these events. 
The second effect is that a late SOS trigger can confuse 
the timing logic in such a way that the coincidence tim- 
ing clock, which usually starts with the HMS and stops 
with the SOS, starts and stops with SOS ("retiming"). 
Wrongly timed events appear at lower (coincidence block- 
ing) and higher (retiming) TDC channel numbers. 

The coincidence blocking factor was calculated for a 
number of runs taken at different trigger rates using 
HMS-SOS raw (not corrected for pathlcngth) coincidence 
time (TDC's) spectra. We found that the coincidence 
blocking correction, k CO i n , depends nearly linearly on the 
rate of the pretrigger, 5-40 kHz for the case of the elec- 
tron spectrometer, the SOS. The coincidence blocking 
correction was then parameterized in terms of 

kemn = 1 - aN strig , (12) 

where a w 2.218 x 10 -5 and N str i g is the SOS trigger rate 
(in Hz). This correction for our coincidence time window 



of 120 ns was up to 4.5%, with an uncertainty of ~0.1%. 

The HMS gas Cerenkov detector is used for electron 
rejection in the tt~ production case. The effective time 
window is given by the Cerenkov detector ADC "gate" 
and is approximately 100 ns wide. The loss of pions due 
to Cerenkov-detector blocking is due to electrons pass- 
ing through the detector after the first particle (pion), 
but within the effective ADC gate window. In this sit- 
uation, the signal from the electron will be associated 
with the original pion trigger and the pion event will be 
mis-identified as an electron. Such mis-identified pions 
arc eliminated due to analysis cuts, so that the elec- 
tron event effectively blocks the HMS Cerenkov detec- 
tor for good pion events. The number of pions lost due 
to the Cerenkov detector blocking depends only on the 
rate of electrons into the spectrometer and does not de- 
pend strongly on variations in run to run characteristics. 
Therefore, we have used a small sample of measurements 
to determine the size of the correction for given kine- 
matics, and for electron rates in the range of 20-500 kHz 
parameterized it in a functional form: 

r cer = l-1.6xHT 7 x^, (13) 

-L run 

where A^ ec / is the clean electron trigger (ELCLEAN) 
counts, defined by high level cuts on calorimeter and 
Cerenkov detector, and T run is the duration of the run 
(in seconds). Cerenkov-detector blocking effect was on 
the level of ~2% at 100 kHz, and reached up to ~6% 
at 400 kHz, with a systematic error less than 1%. The 
uncertainty in the HMS Cerenkov detector blocking cor- 
rection is largely attributed to the uncertainty in the 
Cerenkov detector timing window. In particular the 
effective Cerenkov detector gate width can be slightly 
larger than the measured ADC gate (salOO ns). While 
the ADC gate is fixed, the Cerenkov detector signal itself 
has some width and the overlap determines an effective 
gate width. 

Computer Dead Time: The computer dead time 
strongly depends on the trigger rate and experimentally 
is directly measured by scalers that record the number of 
triggers {N tr ig) and pretriggers {N pretng ). Since pretrig- 
gers are generated for each particle, and triggers are only 
read out for those events for which the Trigger Supervisor 
is not busy, the computer live time is N tr i g /Np retr ig. 

The computer dead time varied from a few percent at 
low rates to up 30% at trigger rate ~2 kHz. The un- 
certainty in the computer live time measurement is esti- 
mated by the deviation of the measured value from the 
value calculated from the total rate. The resulting un- 
certainty is ~ 0.2%. The electronic dead time was always 
< 1% and often negligible. 

For the E00-108 experiment, the computer and elec- 
tronic live time corrections are applied run-by-run. More 
details of the analysis and corrections can be found in 
Ref. [Hil. 

Other corrections: From a measurement detecting 
positrons in SOS in coincidence with pions in HMS, we 
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found the background originating from ir° production 
and its subsequent decay into two photons and then 
electron-positron pairs, or e + e _ 7 directly, negligible. In 
addition, a small ~ 2% correction was made to the deu- 
terium data to account for a small Final-State Interaction 
effect of the pions traversing the deuterium nucleus [H| . 



B. Model Cross Section and Monte Carlo 
Simulations 

We added the possibility of semi-inclusive pion electro- 
production to the general Hall C Monte Carlo package 
SIMC H, usin g E q- ©■ Thc CTEQ5 next-to-lcading- 
order (NLO) parton distribution functions were used to 
parametrize qi(x,Q 2 ) [90j . and the fragmentation func- 
tion parameterization for D+_^(z, Q 2 ) + D~^(z, Q 2 ), 
with D + (D~) the favored (unfavored) fragmentation 
function, from Binnewies et al. [l8j |. The remaining un- 
knowns are the ratio of D~/D + , the slope b of the Pt 
dependence, and the parameters A and B describing thc 
4> dependence. Both thc D~ /D + ratio [9l[ and the b- 
value 92| are taken from HERMES analysis. The latter 
is chosen for consistency with the comparisons shown in 
our earlier publication [l3[ , but closely coincides with the 
averaged value for all data. We will study the detailed 
Pt-dependence of our data later on in Section VIIF. 

When analyzing our data as a function of Pt, we found 
that the Q 2 -dcpendcncc of the cross sections needed to 
be altered slightly from thc factorized high-energy expec- 
tation [sU to obtain a smooth P t dependence. This is not 
too surprising, as the (low) energies of our semi-inclusive 
pion production measurements are beyond the region 
where the BKK fragmentation functions were shown to 
describe existing data. Hence, we introduced an addi- 
tional <5 2 -dcpcndcnt multiplicative term in the model 
cross section in thc form 



TABLE X: Corrections and systematic uncertainties. 



F{Q 2 ) = 1 + C x ■ HQ 2 



9i 

r>2 



9l 
Q 4 ' 



(14) 



The parameters Cx, C2 and C3 were adjusted in such a 
way that the calculated yields from the SIMC simula- 
tion match the experimental data. To accomplish this, 
the ratio of experimental and SIMC yields were calcu- 
lated in a number of Q 2 bins, and the model cross sec- 
tion was iterated until the ratios approach unity. A vari- 
ety of fits, with different combinations of data included, 
more complicated fit functions (including ^-dependent 
terms with additional binning in <p and Pt) rendered 
parameters Cj that remained reasonably stable, within 
±10 — 20%. As average "best values" for the fit parame- 
ters, we adopted d = 0.889, C 2 = -2.902 and C 3 =3.050. 
Recall that for most of the cross section results (at P t « 
0.05 GeV/c) we neglected the ^-dependence and kept the 
parameters A and B at 0, in accordance with both theo- 
retical expectations (discussed in subsection IIIIBI ) , and 
our own findings (see subsection IVII Al ) . 



Source of correction 


Range (%) 


Systematics (%) 


Detector inefficiencies 


5-10 


1-2 


Target wall contribution 


2-3 


1.0 


Accidentals 


10-20 


1-2 


Pion absorption 


1-2 


1.0 


Pion decay 


2-10 


1.0 


Kaon contamination 


0.2-2.0 


0.5 (2 < 0.7) 


Radiative corrections 


5-10 


1-2 


Exclusive tail 


5-15 


0.5-2.5 {z < 0.8) 


Pions from diffractive p 


5-15 


0.5-2.5 


Computer Dead Time 


5-25 


0.2 


Coincidence blocking 


1-4.5 


0.1 


Cerenkov detector blocking 


2-4 


<1 


Other corrections 


1-2 


1.0 


Total 


15-40 


3.5-7.5 



VI. SYSTEMATIC UNCERTAINTIES 

As part of the analysis, several systematic studies were 
performed on the data to verify that the measured cross 
sections and ratios are not biased by the detector, event 
selection and background correction effects. The level of 
corrections applied to the experimental data and related 
systematic uncertainties are listed in Table |X] 

For absolute cross sections we have added all system- 
atic uncertainties in quadrature. Note, that in practice 
the range of applied corrections and related systematic 
errors are slightly different for 7r + and tt~ . For exam- 
ple, the Cerenkov blocking is clearly far larger for n~ . 
Part of the corrections (such as radiative, pion decay, de- 
tector inefficiencies) are nearly identical for ir + and tt~ 
and cancel in the ratios, hence related systematic uncer- 
tainties are much smaller for the ratios. Below we will 
discuss the most dominant sources of systematic uncer- 
tainties related with pions from the radiative tail from 
exclusive pion electroproduction and pions from thc de- 
cay of diffractive p° mesons. 



A. Uncertainties related to the exclusive pion tail 

The model used in SIMC for exclusive pion electropro- 
duction mainly focused on parallel kinematics (with the 
outgoing pion along the direction of the virtual photon) 
for the purpose of understanding the z-scan data. To 
estimate the possible systematic error that arose from 
ignoring the LT and TT interference terms, and to test 
the absolute magnitude of the correction, we extracted 
simulated yields for the exclusive radiative tail calculated 
using our nominal, empirical parameterization as well as 
the MAID model. 

Figures [5l [6] and [7] show the results of these simula- 
tions (MAID = red circles, SIMC= blue squares). Thc 
yields from each model [top panel] as well as the ratio 
between the two [bottom panel] are plotted versus 9 pq . 
While the two calculations differ in the absolute mag- 
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FIG. 5: (Color online) The exclusive radiative tail simula- 
tions for positive pions from deuterium using the nominal 
model cross section in SIMC described in the text. The top 
plot shows the total simulated yield from exclusive radiative 
events vs. the angle between the outgoing pion and the virtual 
photon. All points are at a fixed z = 0.55. For comparison, 
the yield using the MAID model (which has limited validity 
for W > 2 GeV) is also shown. The bottom panel shows the 
ratio of exclusive yields from both models. 



nitudc of the exclusive yield, the ratio of yields shows 
little dependence on the outgoing pion angle. This sug- 
gests that the choice to ignore the interference terms in 
the SIMC parameterization had minimal effect (the con- 
tribution from the interference terms should increase at 
larger pion angles) over the region studied. 

The effective W at the vertex for events that contribute 
to these yields is around 1.9 GeV (on average) so the em- 
pirical parameterization in SIMC, which agrees well with 
the JLab data at W = 1.95 GeV, should be more ap- 
propriate. Half of the difference between the two results 
was used as an estimate of the systematic uncertainty 
for the contribution of radiative exclusive events to the 
semi-inclusive yield. 



FIG. 6: (Color online) The exclusive radiative tail simulations 
for negative pions from deuterium (as in Fig. [5]). 
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B. Uncertainties related to events from diffractive 
p production 



FIG. 7: (Color online) The exclusive radiative tail simulations 
for positive pions from hydrogen (as in Fig. [5]). 



This uncertainty is related to the choice of the param- 
eterization for the p° cross sections. As mentioned, we 
used cross section based on the PYTHIA 82[ generator 



with modifications as implemented by the HERMES col- 
laboration [sH and additional modifications to improve 
agreement with CLAS data 84|. To estimate systematic 



uncertainties related to the diffractive p subtraction, all 



calculations were repeated with slightly (~10%) differ- 
ent values of the parameters. Thus it was found that the 
diffractive p subtraction contributes a systematic uncer- 
tainty of up to ps 2.5% . 
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C. Other sources of systematic uncertainties 

The systematic uncertainties due to normalization, 
target thickness, computer and electronic dead time, 
beam charge measurement, beam energy and spectrom- 
eter kinematics combine to approximately 2%. We note 
that targets and spectrometer polarity were exchanged 
frequently in this experiment, without noticeable effects. 

The overall systematic uncertainty due to spectrometer 
acceptances is estimated to be < 1%. This is because 
the spectrometers have a sufficient wide vertex length 
acceptance to view a 4 cm extended target with limited 
acceptance losses, given that the SOS spectrometer angle 
was limited to relatively forward angles, and that the 
particles of interest cover a central region of the SOS 
momentum acceptance only. The HMS spectrometer has 
a 10 cm uniform vertex length acceptance, and was used 
in the E00-108 experiment to view a 4 cm extended target 
at angles of 20 degrees or less. Hence, the HMS has full 
acceptance. 

It was verified that the level of changes in our results 
due to variation of the values of PID cuts in the analysis 
slightly varies from case to case but is less than 1-2%. 
These variations are within the systematic uncertainties 
assigned to the detector efficiencies, and so are not taken 
as a separate additional source of uncertainty. 

VII. EXPERIMENTAL RESULTS 

Some of the results have been described in two previ- 
ous papers [l3|, [94[. In the first, we observed for the first 
time the quark-hadron duality phenomenon in pion elec- 
troproduction, and the relation with a precocious low- 
energy factorization approach in semi-inclusive deep in- 
elastic scattering. We quantified the latter by construct- 
ing several ratios of pion electroproduction cross sections 
off proton and deuteron targets, and found the ratio of 
favored to unfavored fragmentation functions to closely 
resemble that of high-energy reactions, up to about z = 
0.7 or missing mass M 2 > 2.5 GeV 2 or so. In the sec- 
ond, we studied the transverse momentum dependence 
of semi-inclusive pion production, and found the depen- 
dence from the deuteron target to be slightly shallower 
than from the proton. In the context of a simple model, 
we related these measurements to the initial transverse 
momentum widths of down and up quarks, and the trans- 
verse momentum widths of favored and unfavored frag- 
mentation functions. The results presented here super- 
sede those of Ref . [§4| , which were for reduced statistics 
and improper application of some corrections. 

In this article, we will present our full results in terms 
of cross sections and ratios for semi-inclusive charged- 
pion electroproduction off proton and deuteron targets, 
and relate the findings to the Quark-Parton Model ex- 
pectations, further highlighting the onset of a precocious 
high-energy factorized parton model description. How- 
ever, for completeness we will start with a short Section 



recapitulating two relevant findings of the previous pub- 
lications. 



A. The azimuthal angle (f> and z dependence of the 
cross sections 

In the E00-108 experiment, the average value of the 
angle <f> is correlated with P t , see Fig. @] (the effect is 
tabulated in [sH). However, we will initially show results 
that have only small average value Pt w 0.05 GeV/c, to 
later on come back to the results function of P t , 
including those for P t up to 0.45 GeV/c, in a separate 
Section. 

For the results at low Pt, one expects small to negli- 
gible contribution from the interference terms A and B 
in Eq. [TJ We found no statistically significant difference 
between the results for 7r + or 7r _ , or proton or deuteron 
targets (95|, and therefore combined all four cases to- 
gether. Taking the systematic uncertainties of approxi- 
mately 0.03 into account, we find values of A and B close 
to zero, with no noticeable x or z dependence. When av- 
eraged over all data, we find A = 0.02 ± 0.02 and B = 
-0.04 ±0.02 at P t « 0.05 GeV/c. Folding this back into 
Eq. [TJ we can neglect the azimuthal-angle dependent cor- 
rections to the cross section and ratio results presented 
in the next Sections. The small values of A and B at 
small P t for SIDIS kinematics are consistent with the ex- 
pectations from kinematic shifts due to parton motion as 
described by Cahn J58| and Levelt- Mulders [59j . 

Since the bulk of our data is taken at low P t = 0.05 
GeV/c, we will neglect any ^-dependence and assume A 
= B = until we explicitly revisit the Pt dependence of 
the measured cross sections and ratios later on. 

In our first publication [Hi ] we compared the measured 
1 ' 2 H(c,c , 7r ± )A cross sections as a function of z (at x = 
0.32) with the results of a parton model calculation as- 
suming CTEQ5M Parton Distribution Functions (PDFs) 
at Next-to-Leading Order [9(| and the parameterized 
Fragmentation Functions (FFs) of Binnewies, Kniehl and 
Kramer (BKK) [l8j]. The ratio of unfavored to favored 
fragmentation functions D~ /D + , and the slope b- values 
of the Pt dependences of the cross sections were taken 
from HERMES analysis [H [92|. We found excellent 
agreement between data and Monte Carlo for z < 0.65, 
but striking deviations around z = 0.8. Within our kine- 
matics (at P t ~ 0), M 2 is almost directly related to z, as 
Mj = Mp + Q 2 (l/x - 1)(1 - z). Hence, we attributed 
the large "rise" in the data with respect to the simu- 
lation at z > 0.8 to the N - A(1232) region. Indeed, 
if one considers a 1 H(c,c , 7r _ )A spectrum as function of 
missing mass of the residual system X , one sees only one 
prominent resonance region, the N — A region. Appar- 
ently, above M 2 ps 2.5 GeV 2 , there are already sufficient 
resonances to render a spectrum mimicking the smooth 
z-dependence as expected according to the factorization 
ansatz of Eq. dJJ . 

Much of the data shown later on as function of x, Q 2 
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and Pt, will be centered around z = 0.55, so well within 
the region where effects due to the N-A transition can be 
neglected, and excellent agreement was found between 
the low-energy pion electroproduction data and high- 
energy parton model expectations. We do note that this 
z < 0.7, or M% > 2.5 GeV 2 cut, corresponds to the pre- 
diction of Close and Isgur where duality and low-energy 
factorization would become valid [12, u7f . 



B. The x- and Q 2 -dependence of the cross sections 

At a fixed value z = 0.55, well within the range where 
we found excellent agreement between our cross section 
data and a naive high-energy ansatz in terms of next-to- 
leading-order (NLO) parton distributions (PDFs) con- 
voluted with the BKK fragmentation functions (FFs), 
we will now study the x- and Q 2 -dcpcndcnce of the 
1,2 H(e, e , 7r ± )X cross sections. 

We studied the ^-dependence in the range 0.2 < x < 
0.6 by varying the angle of the scattered electron, while 
keeping the beam energy and the virtual photon energy 
(y « 3.9 GeV) fixed. An additional advantage of this 
choice of z and v is that the corresponding outgoing pion 
momentum is larger than 2 GeV/c, well in the region 
where the ir—N cross sections behave smoothly such that 
final-state interactions do not overly complicate interpre- 
tation of the pion yields. Restricting the kinematics to 
such large pion momenta permits to neglect possible dif- 
ferences in 7r + and tt~ rescattering. In a simple Glauber 
calculation we estimated the total pion absorption correc- 
tion due to rescattering to be 2% for a deuterium target, 
and the difference between ir + and ir~ to be less than 1%. 
We apply the 2% deuteron correction for all 2 H data, and 
assume a constant b = 4.66 (GeV/c)" 2 to describe the 
P t dependence, somewhat different from what we will de- 
rive from the specific Pi-dependent measurements later 
on. Even though this does not affect the cross sections 
represented, it does impact the overall agreement with 
the parton-model calculations, where b comes in as an 
overall normalization. The choice b = 4.66 (GeV/c) -2 is 
chosen in these figures for consistency with the compar- 
isons shown in our earlier publication [l3) . We note that 
this choice is also consistent with the HERMES findings. 

We present in Figs. [5] and |H] a selection of differential 
cross sections for the 1 H(e, e , 7r ± )X and 2 H(e, e'n )X re- 
actions, respectively, at low and high u x se t" values of the 
experiment (by this we mean the x value as calculated 
from the central spectrometer kinematics; using the finite 
spectrometer acceptances we present multiple x-bins). 
Since we vary the scattered electron angle, a variation in 
x set (or x) likewise corresponds to a variation of Q 2 set (or 
Q 2 ). For x set = 0.32, Q 2 set = 2.30 (GeV/c) 2 and for x set 
= 0.53, Q 2 et = 3.8 (GeV/c) 2 cross sections are shown 
along with the model calculations. For simplicity, we 
have only considered CTEQ5M parton distributions at 
NLO (90] an d the BKK 18[ fragmentation functions, al- 
lowing for a slightly modified Q 2 dependence. The scope 
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FIG. 8: (Color online) The 1 H(e,e , 7r+)X (left) and 
1 H(e,e / 7r _ )X (right) cross sections at z = 0.55 as a function 
of x, at QL = 2.30 (GeV/c) 2 ) (top) and at Q 2 set = 3.80 
(GeV/c) 2 ) (bottom), respectively. Solid curves are parton 
model calculations. Solid symbols are data after events from 
diffractive p production are subtracted (see text). 



of this work is to judge how well low-energy pion elec- 
troproduction cross sections (and ratios) compare with 
parton model expectations, and comparisons with other 
possibly more sophisticated model calculations is beyond 
this scope. We conclude that the x-dcpendcncc agrees 
reasonably well with the model calculations, but differ- 
ences in the absolute magnitude of the cross section are 
apparent in certain cases. 

Next, we want to compare the Q 2 -dependence of the 
measured cross sections with the parton model expecta- 
tions. However, as described we varied a; by a change in 
the electron scattering angle, which correlates higher x 
with higher Q 2 . A similar correlation exists within the 
finite spectrometer acceptances. Hence, we need to re- 
move this correlation to present data as function of Q 2 
only, at a fixed value of x. We found that x = 0.40 is 
the optimal value to choose, accessible for each of the 
five settings of the x-scan (x set = 0.26, 0.32, 0.39, 0.46 
and 0.53). The ^-dependence of the parton model, de- 
termined from SIMC simulations over the experimental 
acceptance, was used to scale all data within one a;-scan 
setting to x = 0.40. This was accomplished by using the 
ratios of the normalized yields between data and Monte 
Carlo for each Q 2 bin: Y^ xp jY^ ■ The cross sections of 
different £-scan settings were then corrected to x — 0.40 
using these ratios and the corresponding model cross sec- 
tions at x = 0.40. 

The 1 ' 2 H(e, e , 7r ± )X cross sections for all five settings of 
the x-scan are shown versus Q 2 in Fig.[TUl bin-centered to 
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FIG. 9: (Color online) The 2 H(e,eV+)X (left) and 
2 H(e,e'7r _ )X (right) cross sections at z = 0.55 as a function 
of x, at Q 2 et = 2.30 (GeV/c) 2 ) (top) and at Q 2 Et = 3.80 
(GeV/c) 2 ) (bottom), respectively. Solid curves are parton 
model calculations. Solid symbols are data after events from 
diffractive p production are subtracted (see text). 



x=0.40 with this technique. The curves are the parton 
model calculations, and describe the (3 2 -dcpcndence of 
our data remarkably well. Note that the Q 2 -dependence 
is steeper than naively assumed for a swing in Q 2 from 
about 1.5 to 4.0 (GeV/c) 2 , since these are cross sections, 
"bin-centered" to fixed x = 0.40, which induces trivial 
changes in for instance the beam energy and the resulting 
photon flux. Similarly, one can see that the calculated 
cross section drops very fast around Q 2 = 4.0 (GeV/c) 2 , 
which reflects that one has reached the edge of what is 
kinematically possible with our experimental setup. As 
before, the solid symbols are the data after events from 
coherent diffractive p production are subtracted. Such 
corrections were estimated to be < 10% for these cross 
sections, and do not affect the conclusion that the Q 2 - 
dependence of our data surprisingly conforms to the high- 
energy (quark-parton) expectations. 



energy parton model expectations, we now turn our at- 
tention to various ratios constructed from the data, in an 
effort to quantify the agreement with the quark-parton 
model. Especially the ratio of charged 7r + and tt~ semi- 
inclusive ekctropjxichict^ se ctions (or the ratio of 
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FIG. 10: (Color online) The 1,2 H(e,e'7r ± )X cross sections at 
fixed values of x = 0.40 and z = 0.55, as a function of Q 2 . 
The solid curves are the simple quark-parton model calcu- 
lations following a high-energy factorized description. Solid 
symbols are data after events from diffractive p production 
are subtracted (see text). 



their normalized yields, tt + /tt~) is a quantity relatively 
easy to measure accurately with focusing magnetic spec- 
trometers. In contrast to a large-acceptance detector, the 
acceptance, reconstruction, and detection efficiencies for 
positively- and negatively-charged pions are very similar 
in a focusing magnetic spectrometer, allowing for preci- 
sion comparisons. 



C. The z, x and Q 2 dependence of the cross 
sections ratios (n + /n~ and D/H) 

With the semi-inclusive pion electroproduction data at 
our relatively low energies closely resembling the high- 



With the assumptions of factorization, isospin symme- 
try and charge conjugation (and neglecting heavy quarks 
in the valence-quark region), the cross sections (or nor- 
malized yields) of ir production on protons and neutrons 
at fixed Q 2 can be presented as: 
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a/ (x, z) cx iu{x)D+ (z) + d(x)D~ (z) + Au{x)D- (z) + d{x)D+ (z) 
a p % ~ (x, z) cx iu(x)D~ (z) + d(x)D+ (z) + 4u(x)D+ (z) + d(x)D~ (z) 
a n ' K+ (x,z) cx 4d(x)D+(z) +u(x)D~(z) + 4d(x)D~(z) +u(x)D+(z) 
a n ^~{x,z) oc Ad{x)D (z) +u(x)D+(z) +4d(x)D+(z) + u(x)D~(z), 



with D + and D~ the favored and unfavored fragmenta- 
tion functions, respectively. 

The ratio of charged pion production on proton and 
neutron will then be: 



4u(x)+d(x) + (d(x)+4u(x)) 



(4u(x)+d(x))-^ T ^+d(x)+4u(x) 



T + 4d(x)+u(x) + (u{x)+4d{x))-^^ 



(16) 



{4d{x)+u(x))-^fl+u(x)+4d{x) ' 



It is obvious that the fragmentation functions do not 
completely cancel in the tt + /it~ and D/H ratios, and 
some z-dependence remains carried by the term D~ / D + . 
However, in the ratio of charge-combined cross sections, 
such z-dependence will completely cancel. Some of those 
"super-ratios" were presented in our first publication [l]| , 
and showed validity of the factorization assumption up 
to z ~ 0.65: no z-dependence was found. 

The 7r + /7r~ and D/H ratios versus z: Various ra- 
tios of cross sections of positively- and negatively-charged 
pions and proton and deuteron targets are shown as a 
function of z (at x = 0.32) in Figs. QT] and [TJJ Solid 
(open) circles and squares again represent the data after 
(before) events from diffractive p decay are subtracted. 
We also added existing data for the charged-pion pro- 
duction ratios from Cornell (97j , with the solid and open 
triangles representing data at Q 2 = 2.0 (GeV/c) 2 and x 
= 0.24, and Q 2 = 4.0 (GeV/c) 2 and x = 0.50, respec- 
tively. The solid line is again the simple quark-parton 
model calculation. 

The 7r + /ir~ ratio as measured from the proton tar- 
get is larger than those reported by HERMES |98j |. but 
agrees well with the older Cornell data [97|, and is con- 
sistent (but not equal) to the rise in z as expected from 
the quark-parton model calculation up to z ~ 0.6. At 
values of 0.65 < z <0.85, the ratio decreases because the 
7r~A ++ cross section is larger than the 7r + A° one. The 
sharp rise of the ratio at z > 0.85 is due to exclusive tt + 
production. On the other hand, the tt + /ir~ ratio mea- 
sured on the deuteron reproduces the expected rise from 
the quark-parton model calculation very well. The data 
seem to continue the rising trend for z > 0.7, into the re- 
gion where we noticed effects from the N — A transition 
before. 



In our previous article 13( this was explained within 
the SU(6) symmetric quark model, which essentially re- 
moves the effect of resonance transitions on this partic- 
ular ratio, which is inversely proportional to the ratio of 
unfavored to favored fragmentation functions: D~ / ' D + = 
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FIG. 11: (Color online) The ratio -k + /-k~ for proton (top 
panel) and deuteron (bottom panel) targets as a function of 
z, at x = 0.32. Solid (open) circles and squares represent 
the data after (before) events from diffractive p decay are 
subtracted. Solid and open triangles represent existing Cor- 
nell data [13 at Q 2 = 2.0 (GeV/c) 2 and x = 0.24, and Q 2 
— 4.0 (GeV/c) 2 and x = 0.50, respectively. Stars represent 
HERMES data at average values of <Q 2 )=2.5 (GeV/c) 2 , 
(IF 2 >=28.6 GeV 2 , <i/)=16.1 GeV and (x)=0.082. The solid 
line is a naive quark-parton model calculation. 



(4 — r)/(4r — 1), with r the ratio of tt + over tt~ yields 
off a deuteron target. The observed z-dependence of the 
resulting D~ / D + ratio agreed very well with a fit by the 
HERMES collaboration of their data. In [l3| it was also 
observed that the resulting D~ /D + ratio was indepen- 
dent of x, as it should be, and agreed quite well with 
previous HERMES and EMC data. For completeness 
and future use, we have added in this manuscript the 
D~/D + ratios in tabular format in Table IXj and IXIII 
The columns represent the data before and after events 
from diffractive p decay are subtracted. 

When expressed as a nuclear (deuteron over proton) 
ratio of positively-charged tt + yields and negatively- 
charged 7T~ yields (see Fig. 12), the data show a rela- 
tively flat behavior as a function of z up to about z = 0.7, 
where the N — A transition comes in again. These ratios 
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TABLE XI: The ratio of unfavored to favored fragmentation 
function D~ / D + as a function of z, at x—0.32, evaluated 
at leading order of a s (neglecting strange quarks) from the 
deuterium data. 



z 


D~ 


/D + (after p) 


D~ 


/ D + (before p) 





342 





4620±0 


0710 





4620±0 


0710 





370 





4196±0 


0475 





4449±0 


0465 





398 





4838±0 


0453 





5126±0 


0438 





426 





4764±0 


0429 





5087±0 


0411 





454 





4575±0 


0414 





4940±0 


0392 





482 





4425±0 


0413 





4837±0 


0395 





510 





4059±0 


0318 





4530±0 


0306 





538 





3635±0 


0270 





4134±0 


0257 





566 





3699±0 


0266 





4288±0 


0253 





594 





3638±0 


0274 





4280±0 


0267 





622 





3448±0 


0298 





4124±0 


0284 





650 





3157±0 


0289 





3853±0 


0279 





678 





3587±0 


0314 





4376±0 


0307 





706 





3934±0 


0327 





4800±0 


0319 





734 





3137±0 


0273 





3889±0 


0264 





762 





3164±0 


0254 





3911±0 


0254 





790 





2738±0 


0223 





3375±0 


0223 





818 





2625±0 


0198 





3177±0 


0198 





846 





2380±0 


0177 





2808±0 


0168 





874 





2294±0 


0161 





2607±0 


0159 





902 





2423±0 


0243 





2555±0 


0238 





930 





3025±0 


0739 





2944±0 


0724 



TABLE XII: The ratio of unfavored to favored fragmentation 
function D~/D + as a function of x, at 2=0.55, evaluated 
at leading order of ct s (neglecting strange quarks) from the 
deuterium data. 



X 


D~ 


~/D+ (after p) 


D~ 


/ D + (before p) 





213 





5048±0 


0835 





5682±0 


0800 





238 





4272±0 


0458 





4789±0 


0435 





263 





4008±0 


0341 





4472±0 


0322 





287 





3939±0 


0311 





4361±0 


0294 





312 





4049±0 


0289 





4446±0 


0277 





338 





4278±0 


0285 





4660±0 


0274 





363 





3334±0 


0252 





3631±0 


0242 





388 





3690±0 


0263 





3987±0 


0253 





413 





3476±0 


0262 





3732±0 


0249 





438 





3914±0 


0298 





4177±0 


0287 





463 





3907±0 


0320 





4142±0 


0310 





488 





4198±0 


0362 





4420±0 


0349 





513 





4436±0 


0403 





4546±0 


0395 





538 





4202±0 


0454 





4385±0 


0440 





562 





4721±0 


0581 





4890±0 


0572 





588 





3533±0 


0553 





3668±0 


0539 



appear in reasonable agreement with the quark-parton 
model calculations. 

The x- and Q 2 -dependence of the 7T+/7T and 
D/H ratios: Given that the z-dependence of our low- 
energy semi-inclusive pion electroproduction data show 
a smooth behavior up to z — 0.7, in reasonable agree- 
ment with the quark-parton model expectations, we now 
turn to the x and Q 2 dependence of the various ratios. 
In Figs. Q2] and we show the ratio of positively- to 
negatively-charged pions versus x and Q 2 , respectively, 
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FIG. 12: (Color online) The deuteron over proton (D/H) 
yield ratio for 7r + mesons (top panel) and n" mesons (bot- 
tom panel), as a function of z at fixed x — 0.32. Solid (open) 
symbols represent the new data after (before) events from 
diffractive p decay are subtracted. The solid lines are the 
quark-parton model expectations, plotted up to z — 0.7 where 
effects from the N — A transition may enter the ratios. 



for both proton and deuteron targets. As before, solid 
(open) circles and squares are the data after (before) 
corrections are made to subtract pions originating from 
diffractive p decay. The Q 2 dependence of the deuteron 
to proton ratio for p° production was studied by the 
HERMES collaboration [96[. We have also added exist- 
ing data in Fig. Q2] from Cornell [97j ■ Solid and open tri- 
angles represent data at x = 0.24 and Q 2 = 2.0 (GeV/c) 2 , 
and at x = 0.50 and Q 2 = 4.0 (GcV/c) 2 . As can be 
seen, the Cornell data are in good agreement with our 
data. The positively- to negatively-charged pion ratios 
are also in surprisingly good agreement with the quark- 
parton model prediction. 

The Q 2 dependence of the ratios was extracted by us- 
ing the 7r + and tt~ production cross sections at x = 0.26, 
0.32, 0.39, 0.46 and 0.53, and "bin-centering" these cross 
sections, and thus their ratios, as before, to one common 
value of x = 0.4. This was done by using the hadron part 
of the model cross section in SIMC. The correction was 
checked by running SIMC and taking cross section ratios 
for proton targets at the five x-sc&n central settings. The 
size of the applied corrections amounts to ~ 15% maxi- 
mum for the proton target, and is always below 10% for 
the deuteron target. 

The results, again at a value of z = 0.55, are shown 
in Fig. Q3] The Q 2 dependence of these ratios is in very 
good agreement with the quark-parton model expecta- 
tions, indicated by the solid curve. This teaches that 
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FIG. 13: (Color online) Ratio tt + /n~ for proton (top panel) 
and deuteron (bottom panel) targets as a function of x at 
2 = 0.55. Solid (open) circles and squares are our data af- 
ter (before) events from diffractive p decay are subtracted. 
Solid and open triangle symbols are Cornell data [97| at 
Q 2 =2.0 (GeV/c) 2 and x = 0.24, and Q 2 =4.0 (GeV/c) 2 and 
x = 0.50. The solid lines are the quark-parton model expec- 
tation. 



whereas the Q 2 -dependence f the measured pion elec- 
troprodution cross sections is in reasonable, but not ex- 
cellent, agreement with the quark-parton model expecta- 
tions, as shown in Fig. 1101 any spurious or higher-twist- 
related Q 2 -dependence get completely absorbed in ratios 
(or have an origin in the x-dependence of such ratios, 
as there is a strong kincmatical correlation between x 
and Q 2 within the E00-108 experimental setup). This is 
good news for low-energy access to quark-parton model 
physics in semi-inclusive meson electroproduction. 

Lastly, we show in Figs. [Poland [TH] the deuteron over 
proton ratios for ir + and ir~ electroproduction, as a func- 
tion of x and Q 2 , respectively. These nuclear D /H ratios 
are at a common z = 0.55 again, and the solid curves 
shown correspond as before to the quark-parton model 
expectations. 

The conclusions from these nuclear D/H ratios are 
not unexpected. The dependences on x and Q 2 from 
the quark-parton model is remarkably close to the data, 
again confirming for these ratios that higher-twist ef- 
fects are small or nearly cancel in ratios. The absolute 
magnitudes of the ratios slightly differs from the quark- 
parton model estimates, which reflects the similar differ- 
ence noted in Fig. [T2] Most obvious is the nuclear D/H 
ratio for the ir~ electroproduction case, where the data 
are some 10% higher than the calculated quark-parton 
model ratio. The origin of the discrepancy is not yet 
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FIG. 14: (Color online) The ratio n + /n~ for proton and 
deuteron targets as a function of Q 2 for x = 0.4 and z = 
0.55. Solid (open) symbols are the ratios after (before) yields 
from diffractive p decay are subtracted. The solid lines are 
the simple quark-parton model expectations. 



clear, but on the other hand the data are at a relatively 
large x of 0.4, where the parton distributions themselves 
start having noticeable uncertainties. The latter is inves- 
tigated in more detail in the next section, by constructing 
direct ratios of the d and u valence quark ratios from the 
data. 



D. The ratio of d/u valence quarks constructed 
from charged-pion yields 

The cross section for tt^ 1 production on a deuteron at 
fixed Q 2 can be presented, in similar format and under 
identical assumptions, as the sum of the pion cross sec- 
tions of Eq. (fTST): 



& (x, z) oc (4D+(z) + D-(z))(u(x) + d_(x)) + 
(4D-(z) + D+(z))(u(x) + d{x)) 

a d *~{x,z) oc (4D-(z) + D + (z))(u(x)+d_{x)) + 
(AD+{z) + D-(z))(u(x) + d(x)) 



(17) 



The measured cross sections or yields for ^ produc- 
tion on the proton and deuteron can in the quark-parton 
model be directly used to form relations in terms of the 
u v and d v valence quark distributions: 

<J P * + - a/' cx {D+ - D-){Au v - d v ) 

Od* + - <r d n ~ w ( V + ~ <~ ) + - ^ ) ( 18 ) 

oc (D+ -D-)(3u v +3d v ), 
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FIG. 15: (Color online) The ratio deuteron over proton for 
7r + (top panel) and 7r~ (bottom) as a function of x at z — 
0.55. The solid (open) symbols are the ratios after (before) 
yields from diffractive p decay are subtracted. The solid lines 
indicate the simple quark-parton model calculations. 
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FIG. 16: (Color online) The ratio deuteron over proton for 
7r + (top panel) and n~ (bottom) as a function of Q 2 , at x = 
0.40 and z = 0.55. The solid (open) symbols are the ratios 
after (before) yields from diffractive p decay are subtracted. 
The solid lines indicate the simple quark-parton model calcu- 
lations. 



FIG. 17: (Color online) Top panel: The ratio of valence 
quarks d v /u v as a function of x at z=0.55. Solid circles are 
our data from E00-108 experiment (at P t ~ 0) after events 
from p decay are subtracted. Solid and open squares represent 
data from WA-21/25 [gland EMC [H. Solid triangle sym- 
bols are HERMES data [55l] integrated over the 0.2 < z < 0.7 
range. Bottom panel: The ratio of valence quarks d v /u v as 
a function of z at x=0.32. Solid circles are our data from 
E00-108 after events from diffractive p decay are subtracted. 
The shaded bands on both panels reflect the values of and 
uncertainties in this ratio using CTEQ parton distribution 
functions, based on Eg. [T9l 90f[ . 



where u v = u — u, and d v = d — d. Of course, only 
the full parton distribution u (and d) is physical, but at 
intermediate to large x, x > 0.3, sea quark contributions 
are small and it is common to consider valence quark 
distributions only in this region. 

The d v /u v ratio can be directly extracted from a spe- 
cific combination of the measured proton and deuteron 
tt^~ cross sections as follows: 

(Jp + (x, z) - al (x, z) = 4u v (x) - d v {x) 
1 ~ <jf(x,z)-aj-(x,z) ~ 3[u v (x) + d v (x)Y 



from which one finds 

d v /u v = (4 



3iS,)/(3J£j+l) 



(19) 



(20) 



Studying the x and z (and Pt) dependences of R pd and 
d v /u v thus provides an excellent test of the validity of the 
high-energy factorized view of the SIDIS process, and the 
various assumptions made. 

The ratio d v /u v is shown in Fig. I17[ both as a func- 
tion of x at z—0.55 (top panel), and as a function of z 
at £=0.32 (bottom panel). The ratios extracted from 
our SIDIS data are also compared to WA-21/25 data 
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from neutrino and anti-neutrino deep inelastic scatter- 
ing off proton targets (solid squares) [99[ , and to ratios 
extracted from forward hadron production data from the 
European Muon Collaboration (open squares) (53|. The 
shaded bands on both panels represent the values (includ- 
ing their present uncertainties) as calculated from Eq. [19] 
using CTEQ parton distribution functions [9(| ■ 

The experimentally extracted ratios appear somewhat 
low as compared to the quark-parton model expecta- 
tions using the CTEQ parton distributions, but possi- 
bly within uncertainties. For the results of the present 
experiment, one should not only take into account exper- 
imental systematic uncertainties, but also possible biases 
due to various assumptions in low-energy factorization 
and symmetry in fragmentation functions, etc. Nonethe- 
less, the E00-108 data (at P t ss 0) are in good agreement 
with previous extractions of WA21/25 and EMC, with 
vastly different techniques. 

The undershoot as compared to CTEQ parton distri- 
bution function expectations can be further investigated 
by investigating the dependence on z of the measured ra- 
tios at a fixed value of x (= 0.32). If isospin symmetry 
between favored (D + ) and unfavored (D~) fragmenta- 
tion functions of light quarks (u and d) and anti-quarks 
(u and d) breaks down (D* + + B\ ± Df ^ and 

Dl ^ Dl + ^ Dl + ^ Df ), the ratios of Eq. Umay 
contain additional z-dependent factors, related to asym- 
metries between the fragmentation functions. Thus, a 
dependence of the extracted u d v /u v ratio" on z will be 
a good indication for a breakdown of the symmetry as- 
sumptions, or of the factorized formalism. Indeed, one 
can witness in the bottom panel of Fig. [T7] a sharp in- 
crease of the extracted d v /u v ratio at z > 0.7. This is 
likely not surprising as z > 0.7 corresponds in E00-108 
kinematics to missing mass < 2.5 GeV 2 ), where e.g. 
the A- and higher-resonance contributions become dom- 
inant. 

Below z ~ 0.7, the extracted d v /u v ratio is found to 
be reasonably independent of z, within the uncertainties 
of the data. On average, the data is somewhat low as 
compared to the quark-parton model expectations based 
upon CTEQ parton distribution functions, similar as was 
found in the ^-dependence of this ratio. As a reminder, 
the data presented in Fig. [T7] are at an average low P t ~ 
0.05 GeV/c; we will revisit any possible Pt dependence 
of the extracted ratio later on. 

Even though the extracted d v /u v ratios from the E00- 
108 experiment tend to undershoot the expectations 
based upon CTEQ6 parton distributions, the agreement 
with the existing WA21/25 and EMC data is good, and 
possibly points to the applicability of the assumed fac- 
torization and access to the quark-parton model in rela- 
tively low-energy SIDIS data. This is consistent with our 
earlier findings in Ref. (l3j . 



E. Nuclear Al/D ratios 

We have also analyzed the pion production ratio from 
aluminum to deuterium targets, Al/D, by using the data 
from the "dummy" target cells. The nuclear EMC ef- 
fect, the modification of the (inclusive) nuclear structure 
functions as compared to those of the free nucleon, was 
originally a revelation and firmly injected the subject of 
quarks into nuclear physics. In the valence quark region, 
a linear decrease in the nuclear ratio of structure func- 
tions (typically A/D) of about unity at x = 0.3 to a 
maximum depletion of 10-20% around x = 0.7 has been 
found. For medium to heavy nuclei, A > 12, the effect 
can be well described by either an atomic mass number 
A -1 ' 3 or nuclear density p dependence. 

In semi-inclusive pion production, nuclear effects are 
more complicated, because in addition to influencing the 
electron-quark scattering part, they can affect the quark- 
hadron fragmentation process. For the purpose of the 
present discussion, we assume that the nuclear effects on 
parton distributions and fragmentation functions simply 
factorize. This has by no way been based on rigorous 
experimental verification. 

Experimental results on semi-inclusive leptoproduc- 
tion of hadrons from nuclei are usually presented in terms 
of multiplicity ratios between nuclear (A) and deuteron 
(D) targets as a function of z and v: 



Ka — 



1 



1 



dN^ 



N% IS dz ' Ng IS dz 



dz dz 



(21) 



where the latter applies, if one can ignore EMC-type ef- 
fects, true for x ^0.3. Using the factorized assumption 
and neglecting the nuclear EMC effect for now (we will 
return to the ^-dependence in the nuclear ratios later 
on), we first constructed a nuclear attenuation from the 
ratios of the normalized yields, 



1 A 



dN\ ,dN* 
dz dz ~ Y£ 



(22) 



where Y\ and Yjj are the normalized yields of the electro- 
produced pions from aluminum nuclei and deuterium, re- 
spectively. 

In Fig. [18] we present the ratio of the normalized pion 
electroproduction yields, Al/D, for both ir + (solid circles) 
and 7r~ (solid squares) versus z, at fixed x=0.32. The 
general features of the data, a value of R below unity 
and decreasing with z, are similar to what has been ob- 
served in other experiments and which glob ally h ave been 
explained within various models (see, e.g., [lQOl j and ref- 
erences therein). This applies even in the region of z > 
0.7 where for both the deuteron target and the nearly- 
isoscalar aluminum target nucleon resonances come into 
play (which within the symmetric SU(6) quark model 
cancel out). 

The x-dependences of the Al/D cross section ratio for 
both 7T+ and 7r~, at z = 0.55, are shown in Fig. [19] We 
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FIG. 18: (Color online) The ratio of aluminum over deuteron 
for both 7r + (circles) and 7r _ (squares) data as a function of 
z at x — 0.32. The dashed curve is a calculation bas ed on 
the model for hadron formation of Bialas and Chmaj [l0l| . 
with th e ha dron formation time inserted from a string tension 
model [l02l ]. Th e sol id curve is a prediction based on gluon 
radiation theory |l03| . For the latter, we scaled between data 
for both 14 iV and 64 Cu, assuming (1 - R h A ) ~ A 1/3 , and took 
the average value. 



show data both before and after the events from diffrac- 
tive p production are removed to demonstrate a slightly 
larger impact on the n~ data. The dashed line repre- 
sents the A - dep endent EMC effect fit from the SLAC col- 
laboration [104| fit. The parameterization is normalized 
to take into account hadron attenuation effects. Over- 
all, the agreement is quite good, even if our data scatter 
somewhat and the x-dependence of the cross section ratio 
Al/D is weak in this region. This confirms that nuclear 
ratios already behave like a high-energy parton model 
expectation at relatively low energies and W 2 . 

The Al/D cross section ratios versus Q 2 are extracted 
using again our Monte Carlo simulations to "bin-center" 
the data to a fixed x = 0.40. We show the extracted 
ratios for both ir + and ir~ data at z = 0.55 and x = 0.40 
as a function of Q 2 in Fig. [20] The multiplicity ratio for 
Al/D is nearly flat with Q 2 , comparable to what was also 
observed for the D/H ratio in FiglTHl Similarly, only a 
very weak Q 2 dependence f or this ra t io wa s observed by 
the HERMES experiment [l00|, DM GIU. The dashed 
lines represent constant fits to the data, with a best-fit 
value for tt+ (tt") of 0.556 ± 0.011 (0.520 ± 0.011). 

Recently, there has been dis cussi on on a possible flavor 
dependence of the EMC effect [l07l ] . This would result in 
a possible different depletion of up quarks as compared 
to down quarks. Predictions indicate a somewhat larger 
depletion for up quarks than for down quarks, or equiv- 
alently larger tt + attenuation than 7r~ attenuation. We 
find the opposite, although the uncertainties are large 
and many complicated nuclear effects may contribute, 
including effects that do not obey a factorized form. Fur- 
ther study of this requires precision measurement of the 
z- an d ^-depen dences of these ratios, and their differ- 
ences [l08l . fl09t . 
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FIG. 19: (Color online) The ratio of aluminum over deuteron 
for 7r + (top panel) and ir~ (bottom) as a function of x at 
z = 0.55. Solid (open) symbols are our data after (before) 
events from diffractive p production are subtracted. The cros s 
symbols are Al/D data from the SLAC collaboration |104| ] . 
whereas the dashed lines are the results of their A-dependent 
global fit to nuclear EMC effect. The data and parameteri- 
zation are both normalized to take into account any hadron 
attenuation effects. 



F. The P t dependence of the cross sections 



The extracted cross sections as a function of the pion 
transverse momentum squared P 2 are shown in Fig. [5T] 
and listed in Table 

The solid lines are exponential fits. The acceptance- 
averaged values of cos0 range from -0.3 at low P t to 
nearly -1 at high P t , while the average values of cos2</> 
range from 0.03 at Pt < 0.1 to 1 at high P t . 

A recent study jll0| ana lyzed these data in combina- 
tion with the CLAS data |lll| . and concluded that in 
the kinematics similar to the CLAS data, the Hall C 
data could be relatively well described by a Gaussian 
model with average transverse momentum width of 0.24 
(GeV/c) 2 . The good description of the ^ cross sections 
from different targets was argued to indicate that the as- 
sumption of flavor-independent Gaussian width for both 
the transverse widths of quark and fragmentation func- 
tions was reasonable, in the valence- a; region for z=0.55. 

If taken as standalone data, a careful examination of 
Fig. |2"T1 shows that the P t 2 -dcpendcnces for the four cases 
arc similar, but not identical within statistical uncertain- 
ties. For a more quantitative understanding of the pos- 
sible implications, we study the data in the context of a 
simple model in which the P t dependence is described in 
terms of two Gaussian distributions for each case. 
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FIG. 20: (Color online) The ratio aluminum over deuteron for 
7r + and 7t _ as a function of Q 2 at z — 0.55 and x = 0.40. Solid 
(open) symbols are data after (before) events from diffrac- 
tive p production are subtracted. The solid curves represent 
a constant value of 0.62, as expected from the gluon radia- 
tion calculation at z — 0.55. The dashed lines represent con- 
stant fits to the data, with value for ir + (7r _ ) of 0.556±0.011 
(0.520±0.011). 
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FIG. 21: (Color online) The P t 2 dependence of differential 
cross sections per nucleon for n ± production on hydrogen (H) 
and deuterium (D) targets at («)=0.55 and (a;)=0.32. The 
solid lines are exponential fits. The error bars are statistical 
only. 



The probability of producing a pion with a transverse 
momentum P t relative to the virtual photon (q) direc- 
tion is described by a convolution of the quark distribu- 
tion functions and pt-depcndent fragmentation functions 
D + (z,p t ) and D~ (z,pt), where pt is the transverse mo- 
mentum of the pion relative to the quark direction, with 
the condition P t = zkt + pt assumed. 

Following Ref. [13], we assume that the widths of the 
quark and fragmentation functions are Gaussian and that 
the convolution of these distributions combines quadrat- 
ically. The main difference from Ref. 57J is that we al- 
low separate widths for up and down quarks, and sep- 
arate widths for favored and unfavored fragmentation 
functions. The widths of the up and down distributions 
are denoted by fi u and fid, respectively, and the favored 
(unfavored) fragmentation widths are given by u+ ) . 
Following Cahn [58j and more recent studies [53], we as- 
sume that only the fraction z of the quark transverse mo- 
mentum contributes to the pion transverse momentum. 
We assume further that sea quarks are negligible (typical 
global fits show less than 10% contributions at x = 0.3). 
To make the problem tractable, in the (^-dependence we 
take only the leading-order terms in (Pt/Q), which was 
shown to be a reasonable approximation up to moderate 
P t in Ref. [13] • This simple model then gives: 



a; + (P t ) = C[4 Cl (P t )e- b i^ 



■,-bZPf 



(^)(§r) C2 (P t )e- 
)c 3 (P t )e- b ^ p ? + (£)c 4 (Pt)e- fe ^ 2 



af (P t ) = C{4(±) Ci (P t )e- b « 



-bZPf 



aZ(P t ) = C[4(%)(££)c 2 (P t )e- b * p > 



(23) 



where C is an arbitrary normalization factor, and the 
inverse of the total widths for each combination of quark 
flavor and fragmentation function are given by 



feJ = (zV u + ^)- 1 



(24) 



and we assume <Jd = (<J P + a n )/2. The ^-dependence is 
taken into account through the terms: 



cos{4>)))fi 2 u b+ 
2; 



Ci(P t ) = l + c (P t 
c 2 (Pt) = l + co(P t ,(cos(0)))^6 rf 
c 3 (P t ) = l + C o(P t ,(cos(0))) M 2 u &- 
c 4 (P t ) = 1 + c (P t , (cos(0)) ) M 2 6+ 

c o (P t ,(cos(0)))- 4 ^>^ 



(25) 



(cos(<£)). 



We fit the Pj-dependence of the four cross sections of 
Eq. [23]for the four widths (fj, u , fid, and fx~), C, and 
the ratios D~ / D + and d/u, where the fragmentation ra- 
tio is understood to represent the data-averaged value at 
z = 0.55, and the quark distribution ratio is understood 
to represent the average value at x = 0.3. The fit de- 
scribes the data reasonably well (x 2 = 68 for 73 degrees 
of freedom), and finds the ratio d/u = 0.39±0.03, in good 
agreement with the LO GRV98 fit (HJ for valence quarks 
(about 0.40). The fit also gives a reasonable value for the 
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ratio D-/D+ = 0.43±0.01 (a fit to HERMES results [Hj], 
D-/D+ = 1/(1 + z) 2 , predicts 0.42 at z = 0.55). Both 
d/u and D~ / D + arc largely vmcorrelatcd with the other 
fit parameters and their values are largely determined by 
the magnitude of the cross sections. To estimate the ef- 
fect of experimental systematic uncertainties on our fit 
results, we repeated the fits with: no diffractive p sub- 
traction; 30% smaller exclusive radiative tail subtraction; 
relative target thickness changed by 1%; and difference 
in 7r + and 7r _ absorptions changed by 1%. The last three 
changes had a negligible effect compared to statistical er- 
rors. The first change mainly affected p 2 ,, shifting it to 
a more positive value by almost the size of the statisti- 
cal error, as shown in Fig. [22] We found no significant 
change to the fit parameters upon adding to p^ and p? d an 
average nucleon transverse momentum squared of 0.00 1 
(GeV/c) 2 (evaluated using the Paris wave function |ll2j ) 
for the deuteron model. 

Since the data are at fixed z, the main terms that 
distinguish large fragmentation widths from large quark 
widths are the (^-dependent Cj terms. While there is a 
significant inverse correlation between the two most im- 
portant quark and fragmentation widths, (p u and p +1 
respectively), the fit indicates a preference for p u to be 
smaller than p + as shown in Fig. [22k . The fit also indi- 
cates a preference for pd to be smaller than /i_ as shown 
in Fig. I2"!2b . So in both cases, fragmentation widths ap- 
pear to somewhat dominate over quark widths, within 
our simple model. 

The fit parameters indicate a non-zero k t width 
squared for u quarks (pi = 0.07 ± 0.03 (GeV/c) 2 ), but 
a d-quark width squared that is consistent with zero 
(pi = -0.01 ± 0.05 (GeV/c) 2 ), as illustrated in Fig.[SJb. 
We do note that intrinsic transverse momentum width for 
u and d quarks presented here are far different from the 
earlier published [94| . The previous analysis results used 
a limited and not required cut in the reconstructed vertex 
coordinate reducing statistics, and also has improper cor- 
rections for contributions of pions from both the decay of 
diffractive p production and the exclusive radiative tail. 
Still, the difference in the two results calls for a future 
careful measurement over a large of kinematics (Q 2 , Pt 
and cos(4>)). 

The results are consistent with a di-quark model jl!3l | 
in which the d quarks are only found in an axial di-quark, 
while the u quarks are predominantly found in a scalar 
di-quark. We plotted the results with equal axial and 
scalar di-quarks masses (M a and M s ) of 0.6 GeV; picking 
M a < M s results in p 2 d < p^, and visa verse, with the 
average remaining near 0.06 (GeV/c) 2 . 

Using the fit parameters, we find the magnitude of 
the cos(0) term A &t P t = 0.4 GeV/c to be about 
—0.15 ± 0.05 for all four cases. These results are similar 
in sign and mag nitude to those found in the HERMES 
experiment [114j . 

We find that the fragmentation widths p + and p- are 
correlated, as illustrated in Fig. |2"2"B . although the al- 
lowed range is not large, and the central values (p\ = 
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FIG. 22: (Color online) Fit parameters (crosses) and one- 
standard-deviation contours (continuous ellipses) from the 
seven-parameter fit to the data shown in Fig. 1211 a) u quark 
width squared p^ versus favored fragmentation width squared 
p\; b) p 2 , versus p 2 _\ c) p l u versus p? d \ d) p_ vs p\. The 
dashed and dotted contours are for the case of no diffractive 
p subtraction and a 30% reduction in the size of the exclusive 
radiative tail subtraction, respectively. The larg e do t near 
the middle of panel c is from a di-quark model .113]. The 
dashed straight lines in panels c and d indicate p? u = p 2 , and 
p_ = p\, respectively. 



0.18 ± 0.02 (GeV/c) 2 and p?_ = 0.14 ± 0.02 (GeV/c) 2 ) 
are in reasonable agreement with each other and with the 
fiavor-avcraged value of 0.20 (GeV/c) 2 found in Ref. 
While there is a slight tendency for the favored width to 
be larger than the unfavored one, a reasonable fit can be 
obtained setting the widths equal to each other (x 2 = 71 
for 74 d.f, p\ = p?_ = 0.17±0.03 (GeV/c) 2 ). Taking into 
account the systematic uncertainties, the favored and un- 
favored widths arc consistent with each other. 



G. The P t dependence of the ratios 

7r + /7r~ ratios versus P 2 : The ratios of charged pions 
for proton, deuteron and aluminum targets as a function 
of P t 2 at z = 0.55 and x = 0.32 are shown in Fig. [23] 
Solid (open) symbols are our data after (before) events 
from p decay are subtracted. The solid lines represent 
the expectations from the simple quark-parton model. 

The average values of the pion ratios for deuteron and 
aluminum are smaller than that for the proton, but they 
are nearly flat with P 2 for all three targets. 

D/H ratios versus P 2 : The deuteron over proton ra- 
tios for 7r + (top panel) and ir~ (bottom) as a function of 
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FIG. 23: (Color online) The ratio n + /tt~ for proton, deuteron 
and aluminum as a function of P t 2 at z — 0.55 and x = 0.32. 
Solid (open) symbols are our data after (before) events from 
diffractive p decay are subtracted. The solid lines are simple 
quark-parton model expectations. 
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FIG. 24: (Color online) The ratio deuteron over proton for n + 
(top panel) and ir~ (bottom) as a function of P t 2 at z = 0.55 
and x = 0.32. Solid (open) symbols are our data after (before) 
events from diffractive p decay are subtracted. The solid lines 
are our simple quark-parton model expectations. 



P 2 at z = 0.55 and x = 0.32 are shown in Fig. El Solid 
(open) symbols arc our data after (before) events from p 
decay are subtracted. The solid lines shown correspond 
as before to the simple quark-parton model calculation. 
As can be seen, the D/H ratios for tt + are in good agree- 
ment with the quark-parton model prediction. For Tt~ 
the experimental data on average are ~10-15% higher 
relative to the model expectation. 

P 2 dependence of the d v /u v ratios: For fixed x 
= 0.32 (Q 2 = 2.30 (GeV/c) 2 ) and z = 0.55, we show in 
Fig. [25] the extracted ratios of the down to up valence 
quark distributions d v /u v as a function of _P t 2 . The ex- 
tracted ratios shown before, as function of x and z in 
Fig. [TT1 were from the lowest P t 2 bin only. As before, the 
extracted ratios are on average below the quark-parton 
model expectations, here based upon the GRV parton 
distributions but also consistent with the earlier compar- 
isons with CTEQ parton distributions. Given the statis- 
tical precision of the E00-108 data, it can not be ruled 
out that the d v /u v valence quark distribution ratio may 
have a dependence on P t (or intrinsic quark momentum 
k t ). Such a dependence is in principle possible within 
a transverse-momentum dependent framework [94[. It 
has been calculated to be small for up and down spin- 
averaged parton distributions in Lattice QCD, with far 
larger dep ende nces found in spin-dependent parton dis- 
tributions jl!5j . 

P 2 dependence of the Al/D ratios: In Fig. |2"B1 the 
ratio aluminum over deuteron for 7r + and 7r~ as a func- 
tion of P 2 at x = 0.32 and z = 0.55 is shown. Solid 
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FIG. 25: (Color online) The extracted ratio d v /u v as a func- 
tion of at a;=0.32 and z=0.55. The solid circles are the 
E00-108 data after events from diffractive p decay are sub- 
tracted. The dashed band is a quark-parton model expecta- 
tion using CTEQ parton distribution function parameteriza- 
tions 0. 



(open) symbols are data after (before) events from co- 
herent p production are subtracted. The data show that 
the Al/D ratio is redu ced at high z, as was observed by 
HERMES group jl05| . Our data show slight differences 
in attenuation of 7r + and ir~ . The reduction seems to be 
stronger for 7r + . 

In our kinematic range the multiplicity ratio for Al/D 
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FIG. 26: (Color online) The cross section ratio aluminum over 
deuteron for ir + (top) and n~ (bottom) as a function of P t 2 
at x = 0.32 and 2 = 0.55. Solid (open) symbols are data after 
(before) events from diffractive p production are subtracted. 
The dashed lines are constant fits to the data. 



as a function of P 2 is nearly flat. The dashed lines in 
Fig. l26l represent constant fits to the data, with a best-fit 
value for tt+ and tt" of 0.575±0.010 and 0.538±0.014, 
and a x 2 /ndf of 1.23 and 1.16, respectively. Similar flat 
behavior for the region P 2 <0.2 (GeV/c) 2 was observed 
by HERMES group for variety of nuclei. 

For future use, we have presented in this manuscript 
the various ratios versus z, x and P 2 in tabular format 
in Tables IXlVllXIXl 



VIII. SUMMARY AND CONCLUSIONS 

In summary, we have measured semi-inclusive elec- 
troproduction of charged pions (ir ) from both proton 
and deuteron targets, using a 5.479 GeV energy electron 
beam in Hall C at Jefferson Lab. We have observed, 
for the first time, the quark-hadron duality phenomenon 
in pion electroproduction reactions. This has important 
consequences for a viable access to a quark-parton model 
description in semi-inclusive deep inelastic scattering ex- 
periments at relatively low energies. Several ratios con- 
structed from the data exhibit, provided that W 2 > 4.0 



GeV 2 and z < 0.7 (or beyond the A-resonance region in 
missing mass), the features of factorization in a sequen- 
tial electron-quark scattering and a quark-pion fragmen- 
tation process. We find the azimuthal dependence of the 
data to be small, as compared to the typically larger az- 
imuthal dependences found in exclusive pion electropro- 
duction data, but consistent with data from other groups 
and theoretical expectations [57], HU based on a semi- 
inclusive deep inelastic scattering approach. 

Examination of the Pt dependence of the cross sec- 
tion shows a possible flavor dependence of the transverse- 
momentum dependence of the quark distribution and/or 
fragmentation functions. In the context of a simple model 
with only valence quarks and only two fragmentation 
functions, wc find the transverse momentum k t width 
of u quarks to be larger than that for d quarks, for which 
the width is consistent with zero within the statistical 
uncertainties. We find that the transverse momentum 
p t widths of the favored and unfavored fragmentation 
functions are similar to each other, and both larger than 
the two quark widths. This is consistent with theoret- 
ical expectations based on fits to the world data. Wc 
have shown the sensitivity of our results to be small to 
possible corrections due to both radiative events from 
exclusive pion production channels and pions originating 
from diffractive p scattering (and decay). In many cases, 
the corrections are negligible, although they can become 
large at large values of z. We believe our work will pro- 
vide a fruitful basis for future studies of the quark-parton 
model and more sophisticated model calculations at rel- 
atively low energies. 
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TABLE XIII: Experimental differential cross sections per nucleus (in nb/GeV 3 /c 2 /sr) versus (in GeV 2 /c 2 ) for ty + and n~ 
production on Hydrogen and Deuterium targets. Left (right) part of the table represents the values before (after) events from 
diffractive p decay have been subtracted. Error bars are statistical only. 
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TABLE XIV: The R — -k + /ir ratios versus Z for H, D and Al targets. Left (right) part of the table represents the ratios 
before (after) events from diffractive p decay have been subtracted. Error bars are statistical only. 
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TABLE XV: The R = tv + /tt ratios versus X for H, D and Al targets. Left (right) part of the table represents the ratios before 
(after) events from diffractive p decay have been subtracted. Error bars are statistical only. 
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2 


680±0 


189 


1.554±0.079 


1 


671±0 


265 


2 


803±0 


211 


1.598±0 


087 


i 


782±0 


322 


0.528 


2 


131±0 


162 


1.571±0.089 


1 


554±0 


242 


2 


205±0 


177 


1.613±0 


096 


i 


630±0 


284 


0.544 


2 


061±0 


189 


1.535±0.095 


1 


444±0 


233 


2 


129±0 


206 


1.566±0 


102 


i 


483±0 


259 


0.560 


2 


122±0 


210 


1.454±0.099 


1 


352±0 


247 


2 


191±0 


228 


1.479±0 


106 


i 


377±0 


270 


0.576 


2 


112±0 


217 


1.584±0.122 


1 


061±0 


272 


2 


161±0 


230 


1.613±0 


130 


i 


059±0 


298 


0.592 


1 


946±0 


238 


1.867±0.169 


2 


257±0 


642 


1 


953±0 


240 


1.879±0 


171 


2 


242±0 


643 



J 



r 



j 



r 



j 



r 



j 



[1] E. D. Bloom and F. J. Oilman, Phys. Rev. Lett. 25, [2] P. D. B. Collins, An Introduction to Regge Theory 
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TABLE XVI: The R = tt + /it ratios versus P t 2 (in GeV 2 /c 2 ) for H, D and Al targets. Left (right) part of the table represents 
the ratios before (after) events from diffractive p decay have been subtracted. Error bars are statistical only. 



p'i 




R H ±dR 




R D ±dR 




R M ±dR 




R H ±dR 




R D ±dR 









004 


A 


072±0 


059 


1 

1 


601±0 


039 


1 
1 


638±0 


103 


A 


201±0 


069 


1 

1 


679±0 


045 


1 
1 


/ o4ztU. loo 





012 


1 
1 


940±0 


065 


1 

1 


535±0 


041 


1 
1 


695±0 


119 


A 


058±0 


075 


1 

1 


603±0 


047 


1 
1 


770-1-0. 1 Afi 





020 


1 
1 


951±0 


077 


1 

1 


575±0 


050 


1 
1 


372±0 


101 


■') 
A 


071±0 


090 


1 

1 


653±0 


059 


1 
1 


AAA -l-ft 1 OA 





028 


1 
1 


978±0 


093 


1 

1 


464±0 


054 


1 
1 


436±0 


125 


A 


109±0 


110 


1 

1 


516±0 


063 


1 
1 


^OQ-l-ft 1 fiO 





036 


1 
1 


873±0 


099 


1 
1 


564±0 


064 


1 
1 


642±0 


154 


i 
1 


991±0 


117 


1 
1 


645±0 


076 


1 
1 


QO^-l-fl 01 n 
OZOztU.zlU 





044 


1 
1 


982±0 


113 


1 

1 


508±0 


068 


1 
1 


479±0 


154 


9 
A 


126±0 


136 


1 

1 


573±0 


080 


1 
1 







052 


A 


099±0 


131 


1 
1 


682±0 


085 


1 
1 


621±0 


181 


9 
A 


266±0 


160 


1 

1 


793±0 


104 


1 
1 


am — i— r» oka 





060 


■') 
A 


176±0 


157 


1 
1 


460±0 


078 


1 
1 


361±0 


170 


9 
A 


378±0 


196 


1 

1 


529±0 


093 


1 
1 


AOKA-Ci 00(\ 





068 


A 


056±0 


164 


1 

1 


512±0 


089 


1 

1 


859±0 


232 


9 
A 


239±0 


205 


1 

1 


595±0 


108 


2 


141±0.337 


n 
u 




1 
1 




1 fil 


1 

1 


71 n+n 


1 1 1 


1 
1 


10U3ZU 


1 

1UO 


1 
1 




1 

lyo 


1 

1 


oil znu 




1 


220±0.217 





084 


l 


714±0 


146 


1 


512±0 


105 


1 


677±0 


232 


l 


815±0 


175 


1 


605±0 


129 


1 


891±0.331 





092 


2 


176±0 


207 


1 


528±0 


108 


1 


333±0 


199 


2 


379±0 


262 


1 


617±0 


133 


1 


424±0.266 





100 


2 


722±0 


294 


1 


563±0 


115 


1 


459±0 


223 


3 


098±0 


402 


1 


663±0 


142 


1 


599±0.307 





108 


2 


294±0 


256 


1 


595±0 


127 


1 


497±0 


232 


2 


591±0 


349 


1 


723±0 


164 


1 


628±0.320 





116 


1 


884±0 


215 


1 


649±0 


138 


1 


627±0 


283 


1 


896±0 


263 


1 


786±0 


178 


1 


864±0.425 





124 


1 


897±0 


228 


1 


701±0 


157 


1 


108±0 


200 


2 


069±0 


292 


1 


857±0 


206 


1 


138±0.259 





132 


2 


083±0 


279 


1 


525±0 


154 


2 


023±0 


441 


2 


312±0 


369 


1 


631±0 


197 


2 


527±0.769 





140 


1 


813±0 


258 


1 


645±0 


182 


3 


418±0 


843 


1 


976±0 


332 


1 


795±0 


239 


5 


186±2.065 





148 


1 


896±0 


298 


1 


716±0 


225 


1 


913±0 


543 


2 


073±0 


384 


1 


904±0 


307 


2 


467±1.041 





156 


2 


631±0 


553 


1 


808±0 


260 


2 


135±0 


746 


3 


071±0 


803 


1 


989±0 


343 


2 


871±1.527 





164 


3 


509±1 


355 


2 


088±0 


451 


1 


511±0 


480 


5 


313±3 


421 


2 


533±0 


741 


1 


663±0.660 





172 


2 


754±0 


945 


1 


612±0 


368 


1 


271±0 


607 


3 


439±1 


600 


1 


797±0 


520 


1 


412±0.979 





180 


1 


525±0 


462 


2 


205±0 


639 


4 


001±3 


144 


1 


607±0 


556 


2 


593±0 


962 


10.05±21.82 





188 


1 


049±0 


463 


1 


631±0 


551 


2 


233±1 


506 


1 


058±0 


555 


1 


802±0 


756 


2.917±2.849 





196 


2 


727±2 


279 


4 


250±3 


553 


1 


060±1 


050 


3 


235±3 


431 


6 


571±9 


187 


1.075±1.323 



TABLE XVII: The deuteron over proton (Rd/h~ D/H) and aluminum over deuteron (i? J 4;/ D =Al/D) ratios for tv + and 7r~ 
versus z. Left (right) part of the table represents the ratios before (after) events from diffractive p decay have been subtracted. 
Error bars are statistical only. 



z 


R + D/H ±dR 


R D/H ±dR 


R Al/D ±dR 


R Al/D ±dR 


R + D/H ±dR 


R D/H ±dR 


R Al/D ±dR 


R Al/D ±dR 





317 





814±0 


109 


1 


115±0 


192 





770±0 


146 





763±0 


159 





806±0 


112 


1 


123±0 


207 





759±0 


153 





747±0 


168 





352 





733±0 


034 





927±0 


055 





689±0 


052 





691±0 


060 





721±0 


035 





920±0 


059 





670±0 


055 





663±0 


065 





387 





773±0 


025 


1 


024±0 


050 





622±0 


038 





688±0 


049 





763±0 


026 


1 


025±0 


054 





600±0 


040 





659±0 


053 





422 





747±0 


023 





985±0 


043 





636±0 


039 





597±0 


045 





736±0 


024 





981±0 


046 





614±0 


041 





563±0 


049 





457 





748±0 


023 





995±0 


040 





631±0 


039 





579±0 


043 





737±0 


023 





995±0 


043 





610±0 


041 





543±0 


046 





493 





758±0 


024 





973±0 


037 





602±0 


038 





618±0 


041 





744±0 


025 





970±0 


041 





573±0 


041 





575±0 


046 





527 





702±0 


023 





980±0 


034 





635±0 


038 





562±0 


039 





685±0 


024 





977±0 


039 





605±0 


041 





504±0 


044 





562 





742±0 


024 


1 


004±0 


033 





601±0 


035 





576±0 


038 





727±0 


025 


1 


003±0 


038 





569±0 


037 





520±0 


043 





597 





702±0 


020 


1 


039±0 


032 





606±0 


033 





577±0 


035 





683±0 


021 


1 


041±0 


037 





570±0 


036 





517±0 


040 





632 





778±0 


022 


1 


052±0 


033 





579±0 


032 





635±0 


037 





760±0 


024 


1 


060±0 


039 





536±0 


035 





572±0 


043 





668 





795±0 


021 


1 


016±0 


029 





670±0 


033 





582±0 


034 





774±0 


023 


1 


017±0 


036 





628±0 


037 





494±0 


041 





702 





802±0 


021 


1 


063±0 


031 





607±0 


033 





584±0 


034 





781±0 


023 


1 


077±0 


039 





553±0 


037 





493±0 


040 





738 





882±0 


023 


1 


025±0 


030 





580±0 


031 





526±0 


033 





868±0 


025 


1 


030±0 


038 





524±0 


035 





416±0 


040 





772 





996±0 


024 


1 


004±0 


029 





540±0 


027 





612±0 


034 





996±0 


028 


1 


002±0 


037 





472±0 


031 





507±0 


043 





808 


1 


115±0 


026 





898±0 


024 





543±0 


025 





542±0 


032 


1 


130±0 


030 





875±0 


029 





483±0 


028 





428±0 


039 





842 


1 


044±0 


022 





802±0 


020 





430±0 


021 





457±0 


030 


1 


048±0 


024 





768±0 


023 





371±0 


023 





339±0 


036 





877 





933±0 


023 





973±0 


032 





406±0 


026 





448±0 


037 





927±0 


025 





967±0 


039 





345±0 


028 





329±0 


044 





913 





674±0 


067 





000±0 


000 





382±0 


137 





622±0 


117 





627±0 


075 





000±0 


000 





243±0 


165 





525±0 


146 





947 





693±0 


049 





049±0 


012 


1 


971±0 


119 


5 


758±1 


535 





693±0 


049 





049±0 


012 


1 


971±0 


119 


5 


758±1 


535 





983 





330±0 


025 





211±0 


019 


1 


758±0 


130 


1 


748±0 


214 





330±0 


025 





211±0 


019 


1 


758±0 


130 


1 


748±0 


214 



and High Energy Physics, Cambridge University Press, 
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TABLE XVIII: The deuteron over proton (R D/H 
versus x. Left (right) part of the table represents the 
Error bars are statistical only. 



D/H) and aluminum over deuteron (i?Ai/£>=Al/D) ratios for ty + and ty 
ratios before (after) events from diffractive p decay have been subtracted. 



X 




R D/H ±dR 


R Ai/o ±dR 


R Al/D ±dR 


R D/H ±dR 


R D/H ±dR 




R A;/o ±dR 





208 





770±0 


052 


1 


102±0 


085 





727±0 


066 





742±0 


074 





748±0 


056 


1 


119±0 


103 





693±0 


074 





698±0 


086 





224 





873±0 


057 





963±0 


057 





765±0 


063 





595±0 


058 





859±0 


062 





956±0 


066 





737±0 


069 





527±0 


067 





240 





746±0 


040 





962±0 


051 





639±0 


053 





603±0 


054 





727±0 


043 





955±0 


059 





600±0 


057 





536±0 


063 





256 





768±0 


041 


1 


023±0 


048 





771±0 


058 





539±0 


046 





749±0 


044 


1 


026±0 


055 





747±0 


064 





474±0 


052 





272 





744±0 


037 





993±0 


046 





564±0 


048 





595±0 


047 





726±0 


039 





991±0 


052 





524±0 


052 





537±0 


053 





288 





685±0 


034 





958±0 


043 





673±0 


054 





622±0 


046 





666±0 


036 





947±0 


048 





644±0 


059 





570±0 


052 





304 





706±0 


035 





962±0 


042 





696±0 


057 





547±0 


043 





689±0 


037 





953±0 


047 





669±0 


062 





487±0 


048 





320 





719±0 


036 





983±0 


040 





724±0 


059 





533±0 


041 





703±0 


038 





981±0 


045 





701±0 


064 





477±0 


046 





336 





734±0 


035 


1 


065±0 


046 





670±0 


053 





575±0 


042 





719±0 


037 


1 


072±0 


052 





646±0 


057 





525±0 


047 





352 





730±0 


036 





947±0 


040 





632±0 


052 





546±0 


042 





715±0 


037 





939±0 


044 





606±0 


056 





494±0 


047 





368 





687±0 


031 





908±0 


038 





526±0 


048 





560±0 


042 





673±0 


032 





897±0 


042 





553±0 


053 





510±0 


047 





384 





615±0 


027 


1 


018±0 


044 





536±0 


045 





615±0 


043 





600±0 


027 


1 


020±0 


049 





506±0 


048 





575±0 


047 





400 





757±0 


033 





926±0 


040 





422±0 


039 





558±0 


042 





746±0 


034 





917±0 


044 





388±0 


041 





511±0 


046 





416 





655±0 


028 





987±0 


046 





470±0 


040 





579±0 


044 





641±0 


029 





982±0 


050 





439±0 


042 





537±0 


048 





432 





684±0 


029 





989±0 


046 





544±0 


042 





506±0 


043 





672±0 


030 





988±0 


050 





517±0 


045 





459±0 


047 





448 





725±0 


032 


1 


008±0 


051 





576±0 


045 





564±0 


046 





715±0 


033 


1 


003±0 


055 





553±0 


047 





526±0 


050 





464 





711±0 


032 





981±0 


051 





521±0 


043 





564±0 


050 





700±0 


033 





979±0 


055 





496±0 


045 





526±0 


054 





480 





702±0 


033 


1 


052±0 


062 





579±0 


047 





629±0 


056 





692±0 


034 


1 


056±0 


068 





557±0 


049 





597±0 


061 





496 





633±0 


031 


1 


094±0 


069 





608±0 


051 





607±0 


061 





621±0 


032 


1 


102±0 


076 





588±0 


053 





574±0 


065 





512 





599±0 


031 


1 


002±0 


069 





544±0 


055 





490±0 


064 





588±0 


032 





998±0 


074 





521±0 


057 





449±0 


068 





528 





652±0 


039 





958±0 


069 





555±0 


055 





533±0 


069 





641±0 


040 





954±0 


074 





534±0 


058 





498±0 


074 





544 





781±0 


054 


1 


076±0 


090 





640±0 


066 





618±0 


083 





774±0 


055 


1 


081±0 


096 





625±0 


069 





594±0 


088 





560 





731±0 


054 


1 


092±0 


102 





646±0 


076 





538±0 


088 





723±0 


056 


1 


098±0 


108 





632±0 


079 





512±0 


092 





576 





684±0 


053 





945±0 


095 





403±0 


070 





582±0 


109 





677±0 


054 





942±0 


100 





385±0 


072 





561±0 


114 





592 





833±0 


076 





874±0 


105 





618±0 


087 





334±0 


101 





833±0 


076 





871±0 


105 





617±0 


087 





325±0 


102 



TABLE XIX: The deuteron over proton (Rd/h= D/H) and aluminum over deuteron (7?az/d=A1/D) ratios for ir + and tt~ 
versus P t 2 (in GeV 2 /c 2 ). Left (right) part of the table represents the ratios before (after) events from diffractive p decay have 
been subtracted. Error bars are statistical only. 



Pi 


R D/H ±dR 


R D/J/ ±dR 


R Ai/D ±dR 


R 4i/D ±dR 


R D/H ±dR 


R D/H ±dR 


R Al/D ±dR 


R .4i/D ±dR 





004 





724±0 


017 





922±0 


024 





646±0 


027 





623±0 


031 





709±0 


018 





912±0 


027 





618±0 


029 





573±0 


035 





012 





739±0 


018 





975±0 


029 





633±0 


027 





538±0 


030 





723±0 


018 





972±0 


033 





603±0 


029 





477±0 


034 





020 





781±0 


020 





925±0 


035 





602±0 


027 





698±0 


039 





767±0 


021 





915±0 


039 





569±0 


029 





656±0 


044 





028 





711±0 


020 





973±0 


045 





624±0 


031 





588±0 


038 





692±0 


021 





969±0 


051 





589±0 


034 





530±0 


043 





036 





766±0 


024 





942±0 


050 





550±0 


031 





564±0 


041 





750±0 


025 





934±0 


057 





516±0 


034 





500±0 


046 





044 





741±0 


024 





944±0 


056 





560±0 


034 





567±0 


046 





722±0 


026 





935±0 


064 





519±0 


037 





500±0 


052 





052 





694±0 


024 





898±0 


060 





593±0 


038 





596±0 


053 





674±0 


026 





883±0 


068 





565±0 


041 





527±0 


060 





060 





722±0 


027 





998±0 


076 





590±0 


040 





584±0 


055 





702±0 


029 





986±0 


089 





550±0 


044 





516±0 


063 





068 





757±0 


031 


1 


034±0 


086 





673±0 


046 





565±0 


060 





739±0 


033 


1 


040±0 


101 





639±0 


050 





493±0 


068 





076 





782±0 


034 





899±0 


082 





515±0 


043 





679±0 


077 





765±0 


037 





881±0 


095 





466±0 


047 





619±0 


090 





084 





718±0 


035 





827±0 


077 





720±0 


057 





633±0 


075 





694±0 


037 





799±0 


087 





687±0 


063 





565±0 


086 





092 





738±0 


037 





948±0 


100 





614±0 


054 





642±0 


080 





717±0 


039 





928±0 


117 





571±0 


059 





580±0 


092 





100 





705±0 


037 


1 


181±0 


138 





621±0 


058 





651±0 


083 





681±0 


039 


1 


217±0 


172 





579±0 


064 





589±0 


097 





108 





717±0 


042 


1 


054±0 


127 





660±0 


067 





712±0 


095 





692±0 


045 


1 


066±0 


157 





616±0 


075 





651±0 


113 





116 





723±0 


045 





932±0 


112 





610±0 


067 





638±0 


096 





699±0 


048 





919±0 


132 





563±0 


074 





561±0 


114 





124 





762±0 


055 





886±0 


117 





563±0 


075 





832±0 


121 





739±0 


060 





864±0 


137 





507±0 


083 





795±0 


147 





132 





729±0 


060 





978±0 


142 





717±0 


098 





541±0 


104 





704±0 


064 





973±0 


172 





680±0 


110 





442±0 


123 





140 





769±0 


072 





885±0 


135 





688±0 


108 





446±0 


104 





745±0 


078 





862±0 


160 





918±0 


143 





318±0 


124 





148 





714±0 


074 





794±0 


140 





614±0 


114 





550±0 


138 





687±0 


080 





753±0 


163 





560±0 


129 





432±0 


168 





156 





800±0 


093 


1 


164±0 


264 





551±0 


120 





465±0 


143 





783±0 


100 


1 


208±0 


345 





500±0 


132 





345±0 


170 





164 





843±0 


142 


1 


441±0 


589 





918±0 


217 


1 


269±0 


384 





821±0 


159 


1 


758±1 


195 





905±0 


252 


1 


378±0 


561 





172 





695±0 


127 


1 


092±0 


398 





536±0 


192 





680±0 


265 





662±0 


138 


1 


057±0 


525 





458±0 


221 





583±0 


335 





180 





800±0 


175 





553±0 


198 





661±0 


251 





364±0 


272 





781±0 


190 





484±0 


215 





619±0 


281 





160±0 


345 





188 


1 


157±0 


392 





744±0 


328 





817±0 


357 





597±0 


367 


1 


185±0 


472 





696±0 


376 





790±0 


409 





488±0 


453 





196 





850±0 


352 





545±0 


604 





518±0 


415 


2 


079±2 


119 





836±0 


382 





412±0 


698 





466±0 


456 


2 


849±4 


515 



33 



[6] F. E. Close, F. J. Oilman and I. Karliner, Phys. Rev. [45 

D 6, 2533 (1972). 
[7] F. E. Close and F. J. Gilman, Phys. Lett. B38, 541 [46 

(1972) . 

[8] F. E. Close, H. Osborn, A. M. Thomson, Nucl. [47 

Phys. B77, 281 (1974). 
[9] X. Ji, J.-P. Ma, F. Yuan, Phys. Lett. B597, 299 (2004). [■ 

[10] C. E. Carlson, Physics & Instrumentation with 6-12 [49 
GeV Beams, Thomas Jefferson Laboratory Accelerator 
Facility, Ed. S. Dytman, H. Fenker, P. Roos, p. 17 [50 
(1998) [51 

[11] A. Afanasev, C. E. Carlson, and C. Wahlquist, Phys. 
Rev. D 62, 074011 (2000). 

[12] F. E. Close and N. Isgur, Phys. Lett. B509, 81 (2001). 

[13] T. Navasardyan et al, Phys. Rev. Lett. 98, 022001 [52 
(2007). 

[14] F. E. Close and W. Melnitchouk, Phys. Rev. C 68, 
035210 (2003). 

[15] J. T. Dakin, G. J. Feldman, F. Martin, M. L. Perl, and [53 

W. T. Toner, Phys. Rev. Lett. 31, 786 (1973). 
[16] M. Stratmann and W. Vogelsang, Nucl. Phys. B496, 41 

(1997) ; D. de Florian and R. Sassot, Phys. Rev. D 56, [54 
426 (1997). 

[17] A. S. Raskin and T. W. Donnelly, Ann. Phys. [55 

(N.Y.) 191, 78 (1989). 
[18] J. Binnewies, B. A. Kniehl, and G. Kramer, Phys. Rev. 

D 52, 4947 (1995). [56 
[19] E. L. Berger, Nucl. Phys. B85, 61 (1975). 

[20] E. L. Berger, Proc. of the Workshop on Electronuclear [57 
Physics with Internal Targets, Stanford, CA, January 
5-7, 1987. 

[21] P. J. Mulders, in: R. G. Milner (Ed.), EPIC 2000: 
Proceeding of the 2nd Workshop on Physics with an 
Electron Polarized Light Ion Collider, Cambridge, MA, [58 
2000; |arXiv:hep-ph/0010199V l, 2000. 

[22] W. Melnitchouk, R. Ent and C. E. Keppel, Phys. [59 
Rep. 406, 127 (2005). [60 

[23] W. Melnitchouk, AIP Conf. Proc. 588, 267 (2001). 

[24] J. L. Siegrist et al, Phys. Rev. D 26, 969 (1982). [61 

[25] G. Hanson et al, Phys. Rev. D 26, 991 (1982). 

[26] G. Drews et al, Phys. Rev. Lett. 41, 1433 (1978). [62 

[27] F. E. Close and W. Melnitchouk, Phys. Rev. C 79, [63 
055202 (2009). 

[28] N. Isgur, S. Jeschonnek, W. Melnitchouk and [64 

J. W. Van Orden, Phys. Rev. D 64, 054005 (2001). 
[29] I. Niculescu, Ph.D. thesis, Hampton University, 1999. 
[30] I. Niculescu et al., Phys. Rev. Lett. 85, 1182 (2000). [65 
[31] I. Niculescu et al, Phys. Rev. Lett. 85, 1186 (2000). 
[32] M. Gliick, E. Reya, and A. Vogt, Eur. Phys. J. C 5, 461 [66 

(1998) . 

[33] C. E. Carlson, N. C. Mukhopadhyay, Phys. Rev. D 58, [67 
094029 (1998). 

[34] C. E. Carlson, F. J. Gilman, Phys. Rev. D 7, 2258 [i 

(1973) . [69 
[35] G. Baum et al., Phys. Rev. Lett. 45, 2000 (1980). 

[36] K. Abe et al, Phys. Rev. Lett. 78, 815 (1997). [70 
[37] K. Abe et al, Phys. Rev. D 58, 112003 (1998). 
[38] A. Airapetian et al, Phys. Rev. Lett. 90, 092002 (2003). [71 
[39] R. Fatemi et al, Phys. Rev. Lett. 91, 222002 (2003). 
[40] P. Solvignon et al, Phys. Rev. Lett. 101, 182502 (2008) 
[41] A. Airapetian et al, Phys. Lett. B442, 484 (1998). [72 
[42] B. Adeva et al, Phys. Rev. D 58, 112001 (1998). 
[43] P. L. Anthony et al, Phys. Lett. B493, 19 (2000). [73 
[44] J. Yun et al, Phys. Rev. C 67, 055204 (2003). [74 



A. Calogeracos, N. Dombey, and G. B. West, Phys. Rev. 
D 51, 6075 (1995). 

P. Eden, P. Hoyer, and A. Khodjamirian, JHEP 0110, 
040 (2001). 

T. Sloan, G. Smadja, and R. Voss, Phys. Rep. 162, 45 
(1988). 

K. Ackerstaff et al, Phys. Rev. Lett. 81, 5519 (1998). 

E. A. Hawker et al., Phys. Rev. Lett. 80 , 3715 (1998); 
J. C. Peng et al, |arXiv:hep-ph/98042"88) /l, 1998. 

C. J. Bebek et al, Phys. Rev. Lett. 34 (1975) 759. 
C. J. Bebek, A. Browman, C. N. Brown, K. M. Hansen, 
R. V. Kline, D. Larson, F. M. Pipkin, S. W. Raither, 
A. Silverman, and L. K. Sisterson, Phys. Rev. Lett. 37, 
1525 (1976). 

C. J. Bebek, C. N. Brown, M. S. Herzlinger, 
S. D. Holmes, C. A. Lichtenstein, F. M. Pipkin, 
S. W. Raither, and L. K. Sisterson, Phys. Rev. D 15, 
3085 (1977). 

J. Ashman et al, Phys. Lett. B206, 364 (1988), 
Nucl. Phys. B328, 1 (1989). M. Arneodo et al, Nucl. 
Phys. B321, 541 (1989). 

P. L. Anthony et al., Phys. Lett. B458, 529 

(1999); B463, 339 (1999); B493, 19 (2000). 

K. Ackerstaff et al., Phys. Lett. B404, 383 

(1997) ; B444, 521 (1998); A. Airapetian et al., Phys. 
Lett. B442, 484 (1998). 

J. Adams et al., Phys. Rev. Lett. 92, 171801 (2004); 
S. S. Adler et al., Phys. Rev. Lett. 91, 241803 (2003). 
M. Anselmino, M. Boglione, U. D'Alesio, A. Kotzinian, 

F. Murgia and A. Prokudin, AIP Conf. Proc. 792, 981 
(2005); Phys. Rev. D 71, 074006 (2005); M. Anselmino, 
M. Boglione, A. Prokudin and C. Turk, Eur. Phys. J. 
A 31, 373 (2007). 

R. N. Cahn, Phys. Lett. B78, 269 (1978); Phys. Rev. 
D 40, 3107 (1989). 

J. Levelt and P. J. Mulders, Phys. Rev. D 49, 96 (1994). 
E. L. Berger, T. Gottschalk, and D. W. Sivers, Phys. 
Rev. D 23, 99 (1981). 

K. Oganessyan, H. R. Avakian, N. Bianchi, P. Nezza, 

Eur. Phys. J. C 5, 681 (1998). 

A. Metz, Phys. Lett. B549, 139 (2002). 

A. Bacchetta, M. Diehl, K. Goeke, A. Metz, P. Mulders 

and M. Schlegel, JHEP 0702, 093 (2007). 

D. Boer and P. J. Mulders, Phys. Rev.D 57, 5780 

(1998) ; Nucl. Phys. B564, 471 (2000); D. Boer and 
R. D. Tangerman, Nucl. Phys. B461, 197 (1996). 

V. Barone, Z. Lu and B-Q. Ma, 
|arXiv:hep-ph/0512145| ^l, 2005. 

L. P. Gamberg, G. R. Goldstein, and M. Schlegel, Phys. 
Rev. D 77, 094016 (2008). 

R. Ent, H. Mkrtchyan and G. Niculescu, TJNAF Ex- 
periment E00-108, 2000. 

H. P. Blok et al, Phys. Rev. C 78, 045202 (2008). 

R. Asaturyan et al, Nucl. Instrum. Methods Phys. Res. 
A 548, 364 (2005). 

L. W. Mo and Y. S. Tsai, Rev. Mod. Phys. 41, 205 
(1969). 

R. Ent, B. W. Filippone, N. C. R. Makins, R. G. Mil- 
ner, T. G. O'Neill and D. A. Wasson, Phys. Rev. C 64, 
054610 (2001). 

D. M. Koltenuk, Ph.D. thesis, University of Pennsylva- 
nia, 1999. 

D. Gaskell, Ph.D. thesis, Oregon State University, 2001. 

I. Akushevich, A. Ilyichev, N. Shumeiko, A. Soroko, 



34 



A. Tolkachev, Comp. Phys. Comm. 104, 201 (1997). 
T. V. Kuchto and N. M. Shumeiko, Nucl. Phys. B219, 
412 (1983). 

I. Akushevich and N. M. Shumeiko, J. Phys. G: Nucl. 
Part. Phys. 20, 513 (1994). 

P. Brauel, T. Canzler et al, Z. Physik C 3, 101 (1979). 
V. Tadevosyan et al, Phys. Rev. C 75, 055205 (2007). 
T. Horn et al, Phys. Rev. Lett. 97, 192001 (2006). 
D. Drechsel, S. S. Kamalov, and L. Tiator, Nucl. 
Phys. A645, 145 (1999). 

I. Akushevich, N. Shumeiko, and A. Soroko, Euro. Phys. 
J. C 10, 681 (1999). 

T. Sjostrand, L. Lonnblad, S. Mrenna, and P. Skands, 
|arXiv:hep-ph/0308153kl, 2003. 



P. Liebing, Ph.D. dissertation, University of Hamburg, 
2004. 

C. Hadjidakis, Ph.D. dissertation, Institut de Physique 

Nucleaire Orsay, 2002. 

H. Avagyan, private communication. 

T. Navasardyan, Ph.D. dissertation, Yerevan Physics 

Institute, 2007. 

H. Mkrtchyan, doctoral dissertation, Yerevan Physics 
Institute, 2008. 

M. M. Sargsyan, private communications (2005). 

Hall C Monte-Carlo package for SIDIS, JLab Hall C 



documents, http://www.jlab.org/Hall-C/ 



H. L. Lai, J. Huston, S. Kuhlmann, J. Morfin, F. Olness, 
J. F. Owens, J. Pumpin and W. K. Tung, Eur. Phys. J. 
C 12, 375 (2000). 

P. Geiger, Ph.D. dissertation, Heidelberg University, 
1998. 

B. Hommez, Ph.D. dissertation, University Gent, 2003. 
H. Mkrtchyan, El. J. Contemp. Phys. (Armenian 
Academy Sciences) 43, 266 (2008); Izv. NAN Armenii, 
Fizika 43, 410 (2008). 

H. Mkrtchyan et al, Phys. Lett. B665, 20 (2008). 

H. Mkrtchyan, El. J. Nat. Sciences, NAS RA 1, 12 



[96; 
[97 



[99 
[100 
[101 
[102 

[103 

[104 
[105 



[106 
[107 

[108; 

[109 

[110 

[111 
[112 
[113 

[114 

[115; 



(2009). 

A.G.S. Osborne, Ph.D. dissertation, Glasgow Univer- 
sity, 2006. 

C. J. Bebek, C. N. Brown, R. V. Kline, F. M. Pipkin, 
S. W. Raither, L. K. Sisterson, A. Browman, K. M. Han- 
son, D. Larson and A. Silverman, Phys. Rev. D 16, 1986 
(1977). 

A. Airapetian et al, Eur. Phys. J. C 21, 599 (2001). 

G. Jones et al, Z. Phys. C 62, 601 (1994) 

A. Airapetian et al, Nucl. Phys. B780, 1 (2007) 

A. Bialas and T. Chmaj, Phys. Lett. B133, 241 (1983). 

K. Galmeister and U. Mosel, Nucl. Phys. A801, 68 

(2008). 

A. Accardi, V. Muccifora and H. J. Pirner, 
|arXiv:nucl-th/0211011| /4, 2003. 
J. Gomez et al, Phys. Rev. D 49, 4348 (1994). 
Y. Van Haarlem, for the HERMES collaboration, Pro- 
ceedings of 9th Workshop on Non-Perturbative Quan- 
tum Chromodynamics, Paris, France, June 4-8, 2007. 
J. J. van Hunen, Ph.D. thesis, Utrecht University, 2000. 
I. C. Cloet, W. Bentz and A. W. Thomas, Phys. Rev. 
Lett. 102, 252301 (2009). 

D. Dutta, D. Gaskell and K. Hafidi, JLab Proposal 
PR12-09-004, 2009. 

D. Dutta, J.-C. Peng, I. C. Cloet and D. Gaskell, Phys. 
Rev. C 83, 042201 (2011). 

P. Schweitzer, T. Teckentrup and A. Metz, Phys. Rev. 
D 81, 094019 (2010). 

M. Osipenko et al, Phys. Rev. D 80, 032004 (2009). 
M. Lacombe et al, Phys. Rev. C 21, 861 (1980). 
R. Jakob, P. J. Mulders, and J. Rodrigues, Nucl. 
Phys. A626, 937 (1997). 

Bino Maiheu, Ph.D. thesis, Universiteit Gent, 2006. 
Ph. Hagler, B. U. Musch, J. W. Negele, and A. Schafer, 
Europhys. Lett. 88, 61001 (2009). 



